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2.1 Elliptic flow parameterization. For particles below 4 GeV/ c f 2 = Pip^e~^'^'^/^^^^'^ . 
Above 4 GeV/c the high V2 is used. The systematic uncertainty on the 
V2 is given in the last column |25 

3.1 Centrality bins deflnitions used in this thesis in d+Au and Au+Au col- 
lisions. The minimum and maximum are applied on the reference multi- 
plicity. See text for the definition of reference multiplicity in d-l-Au and 
Au+Au collisions |57 

3.2 Angles from Gaussian fits to the off-diagonal projections |72 

3.3 Expected speeds of sound in different phases of matter [78l[79] [73 

3.4 Values for normalization factors a and b with systematic errors. ZDC 
triggered Au+Au collisions are divided into centrality bins as if they were 
minimum bias events. The cross section numbers in quotes correspond to 
the cross sections as if events were minimum bias but not the actual cross 
sections for these events 85 
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1.1 Diagram of the RHIC collider and the accelerators that feed the beam into 
RHIC. Diagram from [27] 

1.2 Two cutaway views of the STAR detector. The bottom figure shows the 
configuration as of 2001. Figures are from [30] 

1.3 Top Left: Main and Space Trim Magnet coils. Top Right: Poletip Trim 
Magnet coils. Bottom: Diagram of magnet coils. Photos and diagram are 
from [21] 

1.4 Position of the ZDC and DX magnets in the beam pipe. The path of the 
ions, protons and neutrons are shown. Figure is from [39] 

1.5 Diagram of the Central Trigger Barrel. The CTB cylinder is shown along 
with detailed views of the slats. Diagram is from [3S] 

1.6 Diagram of the arrangement of the BBC scintillator tiles. Diagram is 
from [10] 

1.7 Diagram of the STAR TPC from [Mj 

1.8 Cutaway view of the field cage cylinder from [34j- Distances (Typ) are in 




1.9 Cutaway view of the outer subsector pad and wire planes. All dimensions 
are in mm. Diagram is from [3l] 

1.10 The pad plane of one sector. The inner subsector is to the right and the 
outer subsector is to the left. Diagram is from [53] 

1.11 Energy loss distribution of charged particles in the STAR TPC as a func- 
tion of pt for 0.5 T magnetic field. Plot is from [M] 

2.1 Left: Reaction plane and modified reaction plane V2 results as a function 
of Pt in different centrality windows. Right: Comparison of V2 values 
using different flow measurement techniques as a function of centrality 
with 'standard' signifying the reaction plane results and 'u2{4}' signifying 
the 4-particle cumulant results. Flow measurements are from Au+Au 
collisions at ^/snn = 200GeV. Plots are from reference [35] 
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2.2 Ratio of 4-particle cumulant and reaction plane V2- Solid points are mea- 
sured. The line is a fit to a second order polynomial. The open points are a 
conservative estimate of the ratio for the centrality bins without 4-particle 
cumulant V2 measurement. Elliptic flow values are from |,45j |24 

2.3 Detector efficiency for 0-5% most central Au+Au collisions at ^/sNN = 200 
GeV as a function of px- Errors are statistical |24 

2.4 Analysis plots for minimum bias d+Au at ^/snn = 200 GeV/c. Top 
Left: raw signal in Arj (open red points) and background from mixed 
events (black line). Bottom Left: background subtracted signal in Arj (red 
points), background subtracted signal with additional normalized back- 
ground subtracted (blue points), and residual background obtained from 
ZYAl (black line). Top Right: raw signal in A0 (open red points) and 
background from mixed events (black line). Bottom Right: background 
subtracted signal in A0 (red points) and residual background obtained 
from ZYAl. The pt ranges of the trigger particles is 4 < < 6 GeV/c 
and that of the associated particles is 0.15 < px < 4.0 GeV/c. Errors are 
statistical [27 

2.5 Background subtracted A0 distributions for d+Au collisions at ^/sNN = 
200 GeV/c (red) for 4 < p^"^ < 6 and 0.15 < p^''°'' < 4 GeV/c. Left: 
Fit with two Gaussians with centroids fixed at and vr. Right: Compared 
with pp [2S] (blue). Errors are statistical |28 

2.6 Associated particle spectra in d-l-Au collisions at ^/sNN = 200 GeV/c. 
Left: Pt spectra. Right: Z = spectra. Top: Near-side spectra. 
Bottom: Away-side spectra. Errors are statistical |29 

2.7 Comparison of d-l-Au and pp spectra at ^/snn = 200 GeV/c for 4 < 
p^^*^ < 6 GeV/c. Both the near-side and away-side spectra are shown. 

The pp spectra is from |25]. Errors are statistical |30 

2.8 Jet-like correlation peak widths. Left: RMS widths as a function of asso- 
ciated Pt- Open symbols are for pp [25] and filled symbols are for d-l-Au 
■s/snn = 200 GeV/c. Circles are for the away-side width in 0. Squares 
are for the near-side width in (p. Triangles are for the near-side width in 
rj. Right: RMS width (filled) and Gaussian fit sigma (open) as a function 

of trigger pt- Errors are statistical |31 
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2.9 Two-particle correlation functions for 2.5 < p^"^ < 4.0 GeV/c and 1 < 
p^''"^ < 2.5 GeV/c in 10-20% Au+Au collisions at ^JWn = 200 GeV/c. 
Top: Raw correlation function (red). Background from event mixing 
(black). Background from event mixing with elliptic flow added in pair- 
wise by hand (blue). Bottom: Background subtraced 2-particle correlation 
function. The background is normalized by ZYAl before subtraction. Er- 
rors are statistical. Plots for all of the Au-|-Au centrality bins can be found 

in the appendix |33 

2.10 Background subtracted correlation functions in Au-|-Au collisions at ^/snn = 
200 GeV/c for 2.5 < p^"'^ < 4.0 GeV/c and 1.0 < p^''"' < 2.5 GeV/c. 
Left: Our results with |r/^""°'=| < 1.0. Right: PHENIX results from [18] 
with 1 77'^'^*°'^ I < 0.35. Error bars are statistical errors. Histograms represent 

1 sigma systematic uncertainty on the elliptic flow, except that the top 
most histogram on the PHENIX results is -2 sigma systematic uncertainty. 
Shaped bands on PHENIX results are the normalization uncertainty. . |34 

2.11 Away-side correlation amplitudes in the hump {n/3 < \A(f) — 7r| < 47r/9) 
and center (|A0 — 7r| < tt/9) as a function of centrality, {Npart/'^Y^^ shown 
in blue and red respectively. Our data is shown in filled symbols and has 
been scaled down by a factor of 0.35 to account for the rj acceptance. The 
PHENIX data is shown in open symbols. The bars shown the statistical 
errors and the systematic errors due to flow are shown in shaded bands 

and caps for our data and the PHENIX data, respectively. |35 

2.12 Background subtracted correlation functions in Au+Au collisions at y^JJvTV = 
200 GeV/c for 2.5 < p^"'-^ < 4.0 GeV/c and 1.0 < p^''°^ < 2.5 GeV/c. 
Left: Our results with < 1.0. Right: Our results with |j7^'^*°c| < 
0.35. Error bars are statistical. Histograms represent 1 sigma systematic 
uncertainty on the elliptic flow |36 

2.13 (a) Proton to vr"^ and (b) anti-proton to vr" ratios for Au+Au collisions in 

two centralities and d+Au collisions at ^/snn = 200 GeV/c from ^19j. . |37 

2.14 Ionization energy loss plotted as /o5'io(dE/dx) vs. logio{p). The lines 
represent ±0" (±8% in relative dE/dx) bands for protons, kaon, electrons, 
and pions. 170 stands for Bichsel's prediction for 30% truncated mean for 
dE/dx. Figure is from |38 

2.15 Na^, distribution for 3 < px < 3.5 GeV/c and \ri\ < 0.35 from 0-5% 
Au+Au collsions at ^/snn = 200 GeV/c. Positive particles are plotted as 
A^Utt plus 10, negative particles are plotted as -Nct, minus 10. Curves are 
from flt, red is pion, green is kaon, blue is proton, and black is the sum of 
the three. Errors are statistical. Plots for other bins can be found in the 
appendix |39 
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2.16 A^cTtt distribution for 3 < < 3.5 GeV/c and \rj\ < 0.35 from 0-5% 
Au+Au coUsions at ^/snn = 200 GeV/c. Left: Positive particles. Right: 
Negative particles. Curves are from fit, red is pion, green is kaon, blue is 
proton, and black is the sum of the other three. The vertical lines show 
the position of the cuts. Na^ greater than red line gives 95% purity of 
charged pions. Na.^ less than blue line give 50% purity of proton/anti- 
protons. Errors are statistical. Plots for other bins can be found in the 
appendix |41 

2.17 Fits to the cuts on Na.^ as a function of pt for vr^ with 95% purity. Plot is 
for Au+Au 0-5% central collisions at ys/vw = 200 GeV/c with \ri\ < 0.35. 
Errors are statistical |41 

2.18 Identified trigger particle correlations in Au+Au collisions at ^/SNN = 200 
GeV/c. Each group of 3 plots are (from top to bottom) raw signals, raw 
signal minus mixed events, and raw signal minus mixed events and V2- 
Groups from left to right, top to bottom are centralities 50-80%, 30-50%, 
20-30%, 10-20%, 5-10%, and 0-5%. All plots are for Au+Au colhsions with 

3 < p^"^ < 4 GeV/c and 0.15 < p^''""" < 3 GeV/c. Errors are statistical. Es 

2.19 Same as Figure 3.18 except for 2 < pi^'""^ < 3 GeV/c ^4 

2.20 Identified trigger particle correlations in d+Au collisions at ^/sNN = 200 
GeV/c. All plots are for 3 < p^"^ < 4. Left: 0.15 < p^'""" < 3. Right: 
2 < p^***"^ < 3. Top: Raw signals. Bottom: Raw signals minus normalized 
mixed events. Errors are statistical |45 

2.21 Near side spectra for identified protons/anti-protons (blue), charged pions 
(red) and all charged particles (black) in Au+Au collisions at ^/sNN = 200 
GeV/c. Panels are d+Au minimum bias and Au+Au centrality bins 50- 
80%, 30-50%, 20-30%, 10-20%, 5-10%, and 0-5% from left to right and top 

to bottom. Errors are statistical |46 

2.22 Away side spectra for identified protons/ anti-protons (blue), charged pions 
(red) and all charged particles (black) in Au+Au collisions at ^/snn = 200 
GeV/c. Panels are d+Au minimum bias and Au+Au centrality bins 50- 
80%, 30-50%, 20-30%, 10-20%, 5-10%, and 0-5% from left to right and top 

to bottom. Errors are statistical |47 

3.1 Cartoon representations of the expected results for 3-particle azimuthal 
correlations for different physics mechanisms. Top: Unmodified back-to- 
back di-jets. Middle: Jets deflected by radial flow or path length depen- 
dent energy loss. Bottom: Conical emission |50 
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3.2 Results for f4/f| [IS]. A fit to a constant in the px range of 1 < < 2 
GeV/c is used to parameterize V4 in terms of fg. The black line shows a 
fit to a constant over the entire pt range. The red curves are from blast 
wave fits |54 

3.3 TPC acceptance in for particles of 1 < < 2 GeV/c with ±0.5 Tesla 
magnetic field. Positive particles are shown in black and negative particles 
are shown in red. Plots are for 0-5% Au-|-Au collisions. Left: Magnetic 
field with postive polarity. Right: Magnetic field with negative polarity. 
Errors are statistical. Plots for other centralities and for 3 < < 4 
GeV/c particles can be found in the appendix |56 

3.4 Raw 3-particle correlation function for 0-12% central Au-|-Au collisions. |58 

3.5 Left: Raw 2-particle correlation function is shown in red, J2. Normalized 
background from mixed events with v2 and vA added, aBincF2, is shown in 
black. Normalization is done by 3-particle ZYAM. Background subtracted 
2-particle correlation is shown in the minipanel, J2. Right: Hard-soft 
background term, J2 ^-82- Plots are from 0-12% central Au-|-Au collisions 

at ^/sNN = 200 GeV/c. Errors are statistical |60 

3.6 Left: Soft-soft background term, a?hJlp'^. Right: Flow component of 
the soft-soft background term, a?'hT^'^ °^ . Plots are from 0-12% central 
Au-|-Au colhsions at ^/sNN = 200 GeV/c |61 

3.7 Left: Trigger flow correlation up to order Center: Trigger flow corre- 
lation of order f^- Right: Trigger flow correlation with non-flow on soft- 
soft, Jj^. Plots are from 0-12% central Au-|-Au collisions at a/J/vtv = 200 
GeV/c |62 

3.8 Left: Background subtracted 3-particle correlation, J3. Right: Sum of 
hard-soft and trigger flow background terms, J2 ® -B2 + Jj^- Plots are 
from 0-12% central Au-|-Au collisions at ^/snn = 200 GeV/c |64 

3.9 Background subtracted 3-particle correlations in (from right to left, top 
to bottom) pp, d+Au, Au+Au 50-80%, Au+Au 30-50%, Au+Au 10-30%, 
Au-FAu 0-10% and ZDC triggered Au-hAu 0-12%. collisions at y/sj^ = 

200 GeV/c |66 

3.10 Useful for visualizing the on-diagonal elongation and the tilting of the 
near-away peaks. The black line is used to represent the jet-axis. The 
red line is used to represent the trigger particle. The bumps represent the 
near-side and away-side distribution |67 
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3.11 Left: Cartoon with approximate positions of squares used for the yields. 
Right: Average yield in 0.7x0.7 squares centered on the near-side (0,0), 
the away-side (7r,7r), on-diagonal (tt ± 1.42, vr ± 1.42), and off-diagonal 
(7r±1.42, 7r=Fl.42). Npart is the number of participants. The solid errors are 
statistical and the shaded are systematic. The ZDC triggered 0-12% most 
central Au+Au collions point is shifted to the left for clarity. Collisions 

are at y/s^N = 200 GeV/c [67 

3.12 Regions that are projected for the on-diagonal and off-diagonal projections 

are shown in blue and red respectively [68 

3.13 Projections of strips of full width 0.7 radians on the away-side. On- 
diagonal projections (blue) and off-diagonal projection (red) in centralities 
(from left to right, top to bottom) pp, d+Au, Au+Au 50-80%, Au+Au 
30-50%, Au+Au 10-30%, Au+Au 0-10% and ZDC triggered Au+Au 0- 
12% collisions at y^J/vw = 200 GeV/c. Error bars are statistical errors. 
Systematic errors on the off-diagonal projections are shown in the yellow 
band. Dashed black lines are at zero. The on-diagonal projections with 
systematic errors can be found in the appendix |69 

3.14 Difference between in the on-diagonal projection and the side Gaussians 
from fit to the off-diagonal projection. Projections are in strips of full 
width 0.7 radians on the away-side in centralities (from left to right, top 
to bottom) pp, d+Au, Au+Au 50-80%, Au+Au 30-50%, Au+Au 10-30%, 
Au+Au 0-10% and ZDC triggered Au+Au 0-12% colhsions at a/s^ = 200 
GeV/c. Error bars as statistical from the on-diagonal projection only. 
Systematic errors are shown in the yellow band. Dashed black lines 

zero |71 

3.15 Angles from Gaussian fits to the off-diagonal projections. The statistical 
errors (solid) are the fit errors. Systematic errors are shaded. The blue 
line is a fit of the points to a constant, yielding 1.42 ± 0.02 (fit error using 
quadrature sum of statistical and systematic errors on the points). Npart 
is the number of participants. The ZDC triggered 0-12% most central 
Au+Au point is shifted to the left for clarity. The dashed line is at 7r/2. 

The collisions are at ^Js^n = 200 GeV/c [72 

3.16 Background subtracted 3-particle correlations for ZDC triggered 0-12% 
most central Au+Au collisions at ^/snn = 200 GeV/c. Trigger particle is 
3 < < 4 GeV/c. Associated particle transverse momentum from left 
to right top to bottom is 0.5 < pr < 0.75, 0.75 < px < 1, 1 < Pt < 1-5, 

1.5 < pt < 2, and 2 < pt < 3 GeV/c [tS 
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3.17 Away-side projections of background subtracted 3-particle correlations for 
ZDC triggered 0-12% most central Au+Au collisions at ^/snn = 200 
GeV/c. The on-diagonal projection is in blue and off-diagonal projection 
in red. Yellow bands represent the systematic error on the off-diagonal 
projection. Trigger particle is 3 < < 4 GeV/c. Associated particle 
transverse momentum from left to right top to bottom is 0.5 < Pt < 0.75, 
0.75 < Pt < I, I < Pt < 1-5, 1.5 < pr < 2, and 2 < < 3 GeV/c The 
projections are of strips of full width of 0.7 radians |74 

3.18 On-diagonal projections with the side Gaussians from the off-diagonal pro- 
jection removed. The red curves are the Gaussians that were subtracted. 
Yellow bands represent the systematic error on the off-diagonal projection. 
Trigger particle is 3 < < 4 GeV/c. Associated particle transverse mo- 
mentum from left to right top to bottom is 0.5 < pt < 0.75, 0.75 < pt < 1, 
1 < Pt < 1-5, 1.5 < Pt < 2, and 2 < pt < 3 GeV/c The projections are 
of strips of full width of 0.7 radians. Plots are for 0-12% ZDC triggered 
Au+Au colhsions at ^/sNN = 200 GeV/c [75 

3.19 Angles from fits to off-diagonal projections as a function of associated 
particle pt for ZDC triggered 0-12% most central Au+Au collisions at 
y/SNN = 200 GeV/c. Trigger particle is 3 < < 4 GeV/c. The solid 
errors are statistical and the shaded are systematic. The blue line is a fit 
to a constant, yielding yielding 1.36 + 0.03 (fit error using quadrature sum 

of statistical and systematic errors on the points) |76 

3.20 Centrality dependence of the width {a of Gaussians) from fits to away- 
side projections of the background subtracted 3-particle correlations. Top 
Left: Width of the side Gaussians in the off-diagonal projection. Top 
Right: Width of the central peak in the off-diagonal projection. Bot- 
tom Left: Width of the central peak in the on-diagonal projection after 
subtraction of the side Gaussians from the fit to the off-diagonal projec- 
tion. Bottom Right: Difference in the widths of the central peaks done 
as al^.diagonai " ^^of f -diagonal' The uumbcrs indicate the constant fit re- 
suits. The solid errors are statistical and the shaded are systematic. . . l77 

3.21 Centrality dependence of the width {a of Gaussians) from fits to near-side 
projections of the background subtracted 3-particle correlations. Top Left: 
Off-diagonal peak width. Top Right: On-diagonal peak width. Bottom: 



Difference in the peak widths as ^ cr'^on- diagonal - ^^Iff- diagonal- The numbers 
indicate the constant fit results. The solid errors are statistical and the 
shaded are systematic [78 
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3.22 Width (cr of Gaussians) of away-side peaks from the off-diagonal projec- 
tions of the background subtracted 3-particle correlations as a function of 
associated particle pt- Left: Width of the side Gaussians. Right: Width 
of the central Gaussian. Plots are for 0-12% ZDC triggered Au-|-Au at 
^/sNN = 200 GeV/c with 3 < p^"^ < 4 GeV/c. The numbers indicate 
the constant fit results. The solid errors are statistical and the shaded are 
systematic [79 

3.23 Width (a of Gaussians) of near-side peaks from projections of the back- 
ground subtracted 3-particle correlations as a function of associated par- 
ticle pt- Top Left: Off-diagonal peak width. Top Right: On-diagonal peak 



width. Bottom: Difference in the peak widths as ^ al^. diagonal - (^Iff -diagonal- 
Widths are from Au+Au 0-12% ZDC triggered data at ./sj^ = 200 
GeV/c. The numbers indicate the constant fit results. The solid errors 
are statistical and the shaded are systematic |80 

3.24 Three-particle correlations with trigger particles of 4 < p^*^ < 6 (top) 
and 6 < p^'"^^ < 10 (bottom) with associated particles of 1 < pt < 2 
GeV/c for 0-12% ZDC triggered Au+Au colhsions at = 200 GeV/c. 
Left: Background subtracted 3-particle correlation. Right: On-diagonal 
projection (blue) and off-diagonal (red). Error bars are the statistical 
error. Yellow band is the systematic error on the off-diagonal projection. |81 

3.25 Left: Average yield in 0.7x0.7 squares centered on the near-side (0,0), 
the away-side (vr,7r), on-diagonal (vr ± 1.42, vr ± 1.42), and off-diagonal 
(tt ± 1.42, TT =1= 1.42) as a function of trigger particle px- Right: Angles 
from fits to off-diagonal projections. The solid errors are statistical and 
the shaded are systematic. Plots are for 0-12% ZDC triggered Au+Au 
collisions at ^/sNN = 200 GeV/c [82 

3.26 Background subtracted 3-particle correlations for the uncertainty on nor- 
malization factor a. Left: Upper limit on a where it is obtained from 

2- particle ZYAl. Center: Lower limit on a where it is the default from 

3- particle ZYAM minus the difference between the upper limit and the 
default. Right: Off-diagonal projection from default with systematic un- 
certainty from the uncertainty on a shown in shaded band. From top to 
bottom plots are pp, d+Au, Au+Au 50-80%, Au+Au 30-50%, Au+Au 
10-30%, Au+Au 0-10%, and ZDC triggered Au+Au 0-12% colhsions at 
7^ = 200 GeV/c Im 
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3.27 Background subtracted 3-particle correlations for the uncertainty on nor- 
malization factor b. Left: b from 5% of bins for ZYAM. Center: b from 
15% of the bins for ZYAM. Right: Off-diagonal away-side projection from 
default with systematics uncertainty from the uncertainty on b shown 
in shaded band. From top to bottom plots are pp, d-l-Au, Au+Au 50- 
80%, Au+Au 30-50%, Au+Au 10-30%, Au+Au 0-10%, and ZDC triggered 
Au+Au 0-12% collisions at = 200 GeV/c [s? 

3.28 Background subtracted 3-particle correlations for uncertainty on V2- Left: 
V2 from the 2-particle cumulant measurement. Center: V2 from the 4- 
particle cumulant method. Right: Off-diagonal away-side projection from 
default with systematic uncertainty from the uncertainty on f 2 shown in 
shaded band. From top to bottom plots are Au+Au 50-80%, Au+Au 30- 
50%, Au+Au 10-30%, Au+Au 0-10%, and ZDC triggered Au+Au 0-12% 
collisions at ^/sNN = 200 GeV/c [89 

3.29 Background subtracted 3-particle correlations for uncertainty on the jet- 
flow. Left: Jet-flow from 5 GeV/c particles. Left Center: Jet-flow same 
as trigger particle, default. Right Center: Off-diagonal away-side projec- 
tion from default with systematic uncertainty from the uncertainty on jet 
flow shown in shaded band. Right: No jet-flow subtraction. From top 
to bottom plots are Au+Au 50-80%, Au+Au 30-50%, Au+Au 10-30%, 
Au+Au 0-10%, and ZDC triggered Au+Au 0-12% collisions at ^/sj^ = 

200 GeV/c I9I 

3.30 Background subtracted 3-particle correlations for uncertainty on the trig- 
ger particle V4. Left: Trigger particle V4 decreased by 20%. Center: Trig- 
ger particle V4 increased by 20%. Right: Off-diagonal away-side projection 
from default with systematic uncertainty from the uncertainty on trigger 
particle V4 shown in shaded band. From top to bottom plots are Au+Au 
50-80%, Au+Au 30-50%, Au+Au 10-30%, Au+Au 0-10%, and ZDC trig- 
gered Au+Au 0-12% collisions at y/s]^ = 200 GeV/c |92 

3.31 Background subtracted 3-particle correlations for uncertainty on the soft- 
soft term. Left: Soft-soft term constructed from weighed sum of soft-soft 
terms for individual multiplicities. Center: Default. Right: Off-diagonal 
away-side projection from default with systematic uncertainty from the 
uncertainty on soft-soft shown in shaded band. From top to bottom plots 
are Au+Au 50-80%, Au+Au 30-50%, Au+Au 10-30%, Au+Au 0-10%, and 
ZDC triggered Au+Au 0-12% colhsions at ^/sj^ = 200 GeV/c |94 
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3.32 Background subtracted 3-particle correlations for uncertianity on the hard- 
soft on the hard-soft term from bin width effects in ZDC triggered Au+Au 
0-12%. Top Left: Default. Top Center: With hard-soft from finer mul- 
tiplicity bins. Top Right: Difference. Bottom: On-diagonal, left, and 
off-diagonal projections, right. The points are from the default and the 
histograms are from the finer mulitplicity binned hard-soft term |95 

A.l Detector efficiency for charged pions for 3 centralities of d+Au collisions 
(top left) and 9 centralities of Au+Au collisions as a function of px at 

= 200 GeV [lOS 

A. 2 A^a^r distributions with fits in Au+Au collisions at a/saT/v = 200 GeV/c. 
Curves are from fit with red for pion, blue for proton and green for kaon. 
Rows are the centrality bins 70-80%, 60-70%, 50-60%, 40-50%, 30-40%, 
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ABSTRACT 

Ulery, Jason Glyndwr Ph.D., Purdue University, December, 2007. Two- and Three- 
Particle Jet-Like Correlations. Major Professor: Puqiang Wang. 

We present results of 2-particle jet-hke correlations, with high pt triggers and 
identified tt^, p, and p triggers in d-l-Au and Au-|-Au collisions and 3-particle jet- like 
azimuthal correlations in pp, d+Au, and Au-|-Au coUisions at y^s^=200 GeV. These 
results use data from the STAR TPC during RHIC runs II, III, and IV. 

Modifications in 2-particle correlations are observed in Au-|-Au collisions. These 
modifications are not seen in pp or d-l-Au collisions. This demonstrates that the 
modifications are due to final state nuclear effects. 

High Pt protons, anti-protons and charged pions are identified by the relativistic 
rise of dE/dx in the STAR TPC. Correlations of charged hadrons with high pt p, p, 
and TT^ show no discernible difference. The results post challenges to recombination 
and coalescence models which are otherwise very successful in explaining the large 
baryon/meson ratio and the splitting of the elliptic flow at intermediate pt- 

In central Au-|-Au collisions, the away-side 2-particle correlation is significantly 
broadened and even double humped in selective kinematic ranges. Three-particle 
correlations were employed to identify the underlying physics mechanism(s). Results 
in pp, d-|-Au and peripheral Au-|-Au collisions show dijet structure with away-side 
kr broadening. Results in mid-central and central Au-|-Au collisions are consistent 
with a near-side jet and on the away-side a combination of conical emission and large 
angle gluon radiation and deflected jets. The associated pt independent emission 
angle suggests Mach-cone shock waves being the underlying physics mechanism for 
the conical emission. The emission angle is measured to be 1.39 ± 0.01 (stat.) ±0.04 
(sys.) in ZDC triggered 0-12% Au-^Au data. 
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1. EXPERIMENT 
1.1 Relativistic Heavy Ion Collider 

The data used in this thesis is from collisions that were carried out at the Rela- 
tivistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL). This 
collider can in principle collide any nuclei as long as a suitable ion source is available. 
So far it has collided protons (A=l), deutrons (A=2), copper nuclei (A=63) and gold 
nuclei (A=197). The top energies are a/s/vw" = 200 GeV for Au+Au and ^/sNN = 500 
GeV for pp collisions where ^/sNN is the center of mass energy per nucleon pair. Cur- 
rent experimental data has been taken at ^/sNN = 200 GeV for Au+Au, Cu+Cu, 
d+Au, pp, and polarized pp. Additional lower energy runs have been carried out at 
y/SNN = 130 GeV, for Au+Au collisions, at y/s^N = 62.4 GeV in Au+Au, Cu+Cu 
and pp collisions, and ^/sNN = 22 GeV for Au+Au and Cu+Cu collisions. There has 
also been higher energy pp collisions at ^/snn = 400 GeV. There are future plans to 
run pp collisions at the full energy and to do low energy scans {^/sNN ~ 5 — 50 GeV) 
with Au+Au collisions. 

Figure 11.11 shows the chain of accelerators at RHIC. The RHIC beam for ions 
begins with a negative ion source at the Tandem Van de Graaff. The Tandem Van de 
Graaff accelerates ions to 1 MeV/nucleon. The ions are partially stripped of electrons 
in the Tandem Van de Graaff and again on exiting. A bending magnet is then used 
to make a charge selection on the ions. The ions of a particular charge (+32 for Au) 
are then further accelerated by the Booster Synchrotron to 95 MeV/nucleon. After 
exiting the Booster Synchrotron, the ions are further stripped of charge (+77 for Au) 
and then injected into the Alternating Gradient Synchrotron (AGS). In the AGS, the 
ions are accelerated to 10.8 GeV/nucleon (for Au) and then fully stripped on exit. 
Polarized protons (for use in spin studies) are injected from the Proton Linac at 200 
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Figure 1.1. Diagram of the RHIC collider and the accelerators that feed 
the beam into RHIC. Diagram from 



MeV directly into the Booster Synchrotron. Siberian Snakes [28] are used in the AGS 
and RHIC to preserve the polarization. After leaving the AGS, the beam is transfered 
to RHIC. Particles are injected into RHIC with a common magnetic rigidity value 
Bp=81.1141 Tm, where B is the strenght of the magnetic field and p is the radius of 
a charged particle in circular motion in the field. Bp is also equal to the momentum 
perpenduclar to the field divided by the charge. Ions can then be accelerated to a 
maximum magnetic rigidity of 839.5 Tm which corresponds to a magnetic field of 
3.458 T. This gives the maximum kinetic energy of 100 GeV/nucleon for Au ions. 

There are two beam lines throughout the entire acceleration procedure, including 
two Tandem Van de Graaffs. This allows for two different types of ions, one in each 
beam line. With two different types of ions we can get asymmetric collisions such 
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as deuteron and gold which has been run at RHIC. The RHIC ring has a 3.834 km 
circumference and has a total of 1700 superconducting magnets cooled to < 4.6 K. 
There are six interaction points, four of which have dedicated heavy-ion detectors. 
The detectors are STAR, PHENIX, PHOBOS and BRAHMS [29j. 

1.2 STAR Experiment 

One of the two large experiments at RHIC is the Solenoidal Tracker at RHIC 
(STAR) ^0]. The STAR experiment was designed to have a large acceptance (in- 
cluding full, 27T, azimuthal acceptance), high precision tracking, momentum determi- 
nation, and particle identification at midrapidity {\ri\ ~ 1.5). Figure [L2] shows two 
different cutaway views of the STAR experiments. 

The entire detector is surrounded by a 0.5 T solenodial magnet [31]. The magnet 
allows for momentum determination of charged particles. Near the beam pipe resides 
a three layer Silicon Vertex Tracker (SVT) [33j surrounded by an additional layer of 
Silicon Strip Detector (SSD) [33]. These detectors cover full azimuth and \ri\ < 1 in 
pseudo-rapidity. The SVT and SSD were to enhance the measurement of hadrons 
with a short lifetime by providing inner tracking. The primary subdetector of STAR 
is a large volume Time Projection Chamber (TPC) [31] which provides tracking and 
particle identification though ionization energy loss. The TPC has full azimuthal 
coverage and covers \ri\ < 1.8 in pseudo-rapidity. For particle tracking in the forward 
direction (large |?7|) there are two Forward Time Projection Chambers (FTPC) [35j. 
These are radial drift TPCs covering full azimuth and 2.5 < \ri\ < 4 in pseudo-rapidity. 
There is also a Time of Flight (TOF) detector in STAR that provides extended 
particle identification. Currently it only covers a small region but is in the process of 
being upgraded to a full barrel TOF at mid-rapidity. There are two ElectroMagnetic 
Calorimeters (EMC) at STAR. These allow for the measurement of transverse energy 
and for triggering on high pt photons, electrons, and electromagnetically decaying 
hadrons. The two EMCs are the Barrel EMC (BEMC) ^ and the Endcap EMC 




Figure 1.2. Two cutaway views of the STAR detector. The bottom figure 
shows the configuration as of 2001. Figures are from [30] . 
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(EEMC) |37]. The BEMC surrounds the TPC while the EEMC is in front of the TPC. 
These two EMCs provide full azimuthal coverage for the combined pseudo-rapidity 
coverage of — 1 < r/ < 2. There are also additional detectors that are used for trigger 
input. The Central Trigger Barrel (CTB) [38] triggers on the flux of charged-particles 
in the midrapidity region. The Zero-Degree Calorimeters (ZDCs) |39] are used for 
energy determination of neutral particles in the forward direction. The Beam Beam 
Counters (BBC) [IQ] are used to determine luminosity in pp collisions. 

1.2.1 STAR Magnet 

A large mostly solenodial magnet surrounds the STAR experiment and is used 
for particle momentum reconstruction. The specifications required for the STAR 
magnet come from the physics goals. The field has to be large enough to measure the 
momentum of high-energy electron tracks. However the larger the field the higher 
the low momentum cut off for tracks that can be measured in the STAR TPC. The 
field must also be very homogeneous to reduce distortions on the drift electrons in 
the TPC. The magnet can provide a near uniform magnetic field over the range 
0.25 < \Bz\ < 0.5 Tesla. The magnet consists of three types of coils: Main, Space 
Trim and Poletip Trim. Most of the field of the magnet is produced from the Main 
coils. The Trim coils are used to reduce distortions. The coils are cooled by a liquid 
cooling system. The magnet and positions of the coils are shown in Figure 11.31 

Magnetic field mapping is done for all three field components Br, B^p, and Bz- 
This was done in 36 azimuthal points, 57 axial locations. This was measured with a 
steerable array of Hall probes (H] from CERN (European Organization for Nuclear 
Research) and supplemented by NMR measurements. The reproducibility of the 
absolute field is better than ± 0.5 Gauss. For full magnetic field (0.5 T) the maximum 
radial field value is ± 50 Gauss and the maximum azimuthal component is less than 
3 Gauss. 




Figure 1.3. Top Left: Main and Space Trim Magnet coils. Top Right: 
Poletip Trim Magnet coils. Bottom: Diagram of magnet coils. Photos 
and diagram are from |31j- 
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Momentum determination from the magnet is done with the Lorentz Force Law, 

F^qvxB, (1.1) 

where F is the force on a particle with charge q in magnetic field B. Our magnetic 
field is, to a high precision, entirely along the z-axis so our equation becomes, 

F = qvtB (1.2) 

where vt is the velocity perpendicual to the magnetic field. Since our magnetic field 
is along the z-axis the velocity is in the transverse direction. This force provides the 
centripetal force for circular motion, 

F = ^VT. (1.3) 
r 

The two forces must be equal so, 

PtVt 

qvxB — 

r 

Pt = qBr (1.4) 

This gives the transverse component of the particle momentum. The axial component 
and total momentum can be found using the measured angle of the particle with 
respect to the beam (z-axis), 9. 

Pt 



P = 



sin{9) 



1.2.2 Trigger Detectors 



The purpose of the trigger detectors is to enable event selection criteria to be 
applied at a rate greater than that at which the slow detectors operate. This is done 
because RHIC has a crossing rate of about 10 MHz and the slow detectors can only 
operate at rates of about 100 Hz. The fast detectors must be used to provide a 
rate reduction of 5 orders of magnitude and intelligently select desired events. Three 
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triggering detectors were used for the data contained in this thesis. They are the Zero 
Degree Calorimeters (ZDCs), the Central Trigger Barrel (CTB), and the Beam-Beam 
Counters (BBCs). In Au+Au minimum bias collisions, cuts are made on the signals 
in both ZDCs (east and west) and the CTB. There is also a cut on the primary event 
vertex obtained from the ZDCs. The ZDC cuts required a coincidence between the 
two ZDCs of summed signal greater than about 40% of a single neutron signal. The 
CTB cut is used to reject nonhadronic events (which removes events with a very low 
number of charged particles in the CTB). This cut rejects some of the desired events 
that are very perpherial and is why we do not use centralities below the 80% most 
perpherial. In central Au+Au collisions, there are much higher cuts on the ZDCs and 
the CTBs. There is a cut on the primary vertex obtained from the BBC. The cuts are 
tuned such that the events taken are about 10% of the total cross section and such 
that the multiplicity distribution matches the minimum bias distribution for the top 
5% most central collisions. In d+Au minimum bias collsions, the trigger cut was on 
the east (Au side) ZDC only. For pp collisions there is a cut on the BBC signal. 

Zero Degree Calorimeters 

The Zero Degree Calorimeters detect neutrons and measure their total energy in 
a cone about the beam line to both sides of the collisions. The total energy is used 
to calculate the neutron multiplicity. The neutron multiplicity is correlated with the 
event geometry and can be used for centrality determination. 

The ZDCs are hadron calorimeters. Each one consists of three modules. The 
modules have a series of tungsten plates with layers of wavelength shifting fibers that 
route Cerenkov light [42j to a photomultiplier tube. The Cerenkov light is produced 
by the shower particle moving though the optical fibers. The DX dipole magnets 
bend charged particles (protons and ions) away from the ZDCs; leaving only neutral 
particles (neutrons) in the ZDCs. The ZDCs make up a 2 mrad cone about the beam 
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Figure 1.4. Position of the ZDC and DX magnets in the beam pipe. The 
path of the ions, protons and neutrons are shown. Figure is from |39j . 



direction and are about 18 meters from the interaction point. The setup of the ZDCs 
is shown in Fig. II. 4[ 

Central Trigger Barrel 

The CTB measures the charged particle multiplicity at midrapidity (|?7| < 1). 
The CTB consists of four cylindrical bands with 60 scintillator slats each. Each band 
covers the full azimuth and one-half unit in pseudo-rapidity. The CTB is positioned 
just outside the TPC at about 4 meters from the beam pipe. Each one of the slats 
consists of a scintillator, a light guide and a photomultiplier tube. Two slats are 
arranged end to end in an aluminum tray for mounting and handling. Figure 11.51 
shows the setup of the CTB. 

Beam-Beam Counters 

The STAR BBCs are scintillator tiles mounted outside of the pole tip magnets. 
There are 16 small tiles near the beam pipe surround by 16 larger tiles on each 
side. The tiles are hexagon in shape. The smaller tiles can be inscribed in a circle 
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Figure 1.5. Diagram of the Central Trigger Barrel. The CTB cylinder is 
shown along with detailed views of the slats. Diagram is from [38] . 
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with 9.64 cm diameter. The larger tiles can be inscribed in a circle with 4 times 
the diameter of the smaller tiles. The scintillators are connected to photomultiplier 
tubes. The arrangement of the tiles for each of the BBCs is shown in Figure 11.61 
When the pp data used in this thesis was taken a total of 18 of the smaller scintillator 
tiles were installed and used. These were connected to 8 photomultiplier tubes. The 
multiplicities determined from the BBCs are used for pp luminosity measurements 
and triggering. 



STAR Beam-Beam Counter Schematic 
Front View 




Figure 1.6. Diagram of the arrangement of the BBC scintillator tiles. 
Diagram is from [40j . 



1.2.3 Time Projection Chamber 

The TPC is the primary tracking device in the STAR detector. It provides track- 
ing, momentum measurement, and particle identifications. The momentum mea- 
surements come from the curvature of the tracks in the magnetic field as previously 
discussed. The particle identifications is accomplished though measurements of the 
ionization energy loss (dE/dx). The detector covers full azimuth and \ri\ < 1.8. 
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Figure 1.7. Diagram of the STAR TPC from [34j- 

Figure [TTTl shows a diagram of the TPC. The TPC sits in the middle of the STAR 
magnet. It is 4.2 m long and has inner and outer diameters of 1 and 4 m, respectively. 
The TPC volume is filled with PIO gas (90% argon and 10% methane). The gas is 
maintained to a high purity by the TPC gas system. This is required because water 
and oxygen in the gas will absorb electrons. The oxygen is kept below 100 parts 
per million and the water less then 10 parts per million. At this level of purity the 
absorption of electons is only a few percent. The transverse diffusion (a = \/2Dt 
where D is the diffusion coefficient and t is the drift time) of the gas is 230 /im / ^J cm 
at full field (0.5 T). For an electron drifting 2.1 m this give a transverse drift of 3.3 
mm. The drift velocity of the gas is 5.45 cm/ ^s. There is a longitudional drift spread 
of about 230 ns (FWHM). The diffusion sets the scale for the readout. 
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Figure 1.8. Cutaway view of the field cage cylinder from Distances 
(Typ) are in mm. 



A uniform electric field is applied to the TPC which is what causes the electrons 
to drift to the ends. The electric field is about 135 V/cm. The electric field is 
provided by a thin high voltage conductive Central Membrane at the center of the 
TPC, the inner and outer concentric field-cages cylinders and readout end caps on the 
ends of the TPC. Uniformity of the electric field is crucial to achieve submillimeter 
track reconstruction precision for electron drift paths of up to 2.1 m. The Central 
Membrane (CM) is a disk with a central hole. The CM is a cathode and kept at a 
potential of 28 kV. It is made of 70 /im of carbon-load Kapton film. It is mounted 
inside the outer cage cylinder and the inner cage cylinder runs though the hole in 
the center. Thirty six narrow stripes of aluminum have been attached to each side 
to provide targets for lasers for calibration. The field cages keep the electric field 
uniform and provide containment for the gas. As shown in Figure 11.81 the field cage 
cylinders consist of two layers of metal coated Kapton. The layers are separated 
by a NOMEX honeycomb. The metal is etched into 10 mm strips with a 1.5 mm 
separation so that the required voltage difference between rings can be maintained. 
This design was optimized for reduced mass, minimization of track distortions from 
multiple Coulomb scattering and reduction of background from secondaries. 

The TPC end-caps contain anodes and pad planes. The readout planes are multi- 
wire proportional chambers (MWPCs) with pad readout. The chambers consist of 
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Outer Sector 
715.14 



Figure 1.9. Cutaway view of the outer subsector pad and wire planes. All 
dimensions are in mm. Diagram is from [34j. 

three wire planes and a pad plane. Figure 11.91 shows a cutaway view of the positions 
of the wires and the pads. The amplification and readout are done by the anode wire 
planes which consist of 20 fim wires. The wire direction is chosen to best determine 
the momentum of very high transverse momentum particles whose tracks do not 
curve much in the magnetic field. This places the anode wires roughly perpendicular 
to the radial direction. In the other direction the resolution is limited by the wire 
spacing (4 mm). The pad dimensions are also optimized to the best position resolution 
perpendicular to radial tracks. The width of the pad is chosen such that the induced 
charge from an avalanche point on the anode wires shares most of its signal with 
three pads. The outer radius pad subsectors have continuous pad coverage to achieve 
the best possible energy loss (dE/dx) resolution. The outer radius subsectors are 
arranged on a rectangular grid with a pitch of 6.7 mm along the wires and 20.0 
mm perpendicular to the wires and a 0.5 mm gap between pads. There is a 4 mm 
separation between the pad plane and the anode wires. The inner subsector pads are 
optimized for good two-hit resolution due to the high track density in that region. 
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The pads in the inner subsector are 3.35 mm along the wire and 12 mm perpendicular 
to the wires. The inner and outer pads for a sector is shown in Figure 11.2.31 The 
distance to the anode wires is reduced to 2 mm to put most of the signal on three 
pads as it is in the outer sector. The smaller pads give better tracking of the low 
momentum particles. The smaller pads in the inner sector however require the use 
of separate pad row instead of continuous coverage due to constraints on the front 
end electronics. The anode wires voltage is set independently for the two sectors to 
maintain a 20:1 signal to noise ratio for tracks from the center of the TPC. A ground 
grid plane of 75 wires is used to terminate the field in the avalanche region and 
provide additional shielding for the pads. The anodes, pads, and grounding grid make 
up the MWPC. 



Outer Pads Inner Pads 

6.2 mm X 19 5 mm 2.85 mm x 11 .5 mm 

Total of 3,942 Pads Total of 1 ,750 Pads 




Figure 1.10. The pad plane of one sector. The inner subsector is to the 
right and the outer subsector is to the left. Diagram is from |34j . 

The outermost wire plane is the gating grid and is located 6 mm from the ground 
grid. The gating grid is used to control the entry of electrons from the TPC drift 
volume into the MWPC. It also serves to keep positive ions from the MWPC from 
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entering the TPC drift volume. When the grid is open it is transparent to the drifting 
electrons. The grid opens when the event is being recorded by putting all wires at 
the same voltage. The grid goes to alternating positive and negative voltages (±75 
V) when the grid is closed. Positive ions from the MWPC are too slow to go though 
during the open period and are captured during the closed period. The combined 
delays of the trigger and the opening time of the gating grid reduce the active length 
of the TPC from 210 cm to 198 cm due to electrons lost before the grid is opened. 

The MWPCs and gating grids are on 12 sectors on each end of the TPC. These 
sectors are arranged around as circle. There is 3 mm of dead space between the 
sectors. Because of the dead space there is an azimuthal dependence to the particle 
detection efficiency. 

The X and y coordinates of a cluster are reconstructed assuming a Gaussian dis- 
tribution for most tracks. Tracks with a large crossing angle are reconstructed using 
a weighted mean algorithm because ionization is deposited on many pads. The z 
coordinate of a cluster is determined by measuring the drift time of the secondary 
electrons from their origin to the anodes and dividing by the average drift velocity. 
The arrival time of the cluster is calculated by measuring the arrival time of electron 
in discrete time intervals of about 100 ns each. The charge weighted average time 
interval is used for the arrival time. The drift velocity must be known to high preces- 
sion for accurate z position reconstruction. The drift velocity changes with pressure 
and small changes to the gas composition. This is minimized by setting the cathode 
voltage so the electric field in the TPC is at the peak in the velocity vs. electric 
field / pressure curve. The peak is broad and flat so small pressure changes will have 
little change on the drift velocity. The drift velocity is independently measured every 
few hours using tracks created from laser beams. The coUision time can be offset by 
trigger delay, drift time from the gating grid to the anode wires, and the shaping of 
the signal in the front end electronics. The timing offset can be adjusted by recon- 
structing the interaction vertex using data from each side of the TPC separately and 
matching them. 
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The positions of secondary electrons are distorted by non-uniformities in the elec- 
tric and magnetic fields and from global misalignments. The typical size of the distor- 
tions is <1 mm. This size of distortion can have an effect on the transverse momentum 
determination of high px particles. To correct for the distortions the magnetic field 
was measured using Hall probes and NMR probes. The electric field was calculated 
from the geometry. The hit position distortions are calculated and corrected for. Af- 
ter correction the point to track fit error is about 50 /xm and the absolute error on a 
point is about 500 /im. 



dE/dx Vs. P n 




Figure 1.11. Energy loss distribution of charged particles in the STAR 
TPC as a function of pt for 0.5 T magnetic field. Plot is from |34j . 
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The energy loss in the TPC gas can be used for identifying particles at low pt 
(below about 1.2 GeV/c) and at high px (above about 3 GeV/c). For the region in 
between the dE/dx bands cross so identification from energy loss is not possible. The 
energy loss as a function of pt is shown in Fig. 11.111 A resolution of 7% in relative 
dE/dx is required to distinguish protons and pions up to 1.2 GeV/c. A resolution 
of 8% in the relative dE/dx have been achived in the data. The resolution depends 
on the gas gain which is pressure dependent. The pressure is kept at 2 mbar above 
atmospheric pressure so it is time dependent. A wire chamber with a ^^Fe source 
measurses the gas gain. Local variations are calibrated by averaging at the pad-row 
level. The energy loss is measured on up to 45 padrows for each track. This has 
too few points to average out ionization fiuctations. Beacuse of this a most probable 
energy loss instead of an average energy loss is used for particle identification. 
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2. TWO-PARTICLE JET-LIKE CORRELATIONS 

Jets make a good probe because their properties in elementary collisions can be cal- 
culated reliably by pQCD so they arc well calibrated and in heavy ion collisions they 
interact with the created medium. Since they interact with the medium we can study 
the effect of the medium on jets and the effect of jets on the medium. In heavy-ion 
collisions it is not possible to reconstruct jets event by event due to large backgrounds. 
To study jets in heavy- ion collisions we reconstruct jet-like correlations statistically 
through angular correlations. In two-particle correlations this is done by triggering 
on an intermediate or high pt particle and studying the angular distributions of the 
other particles in the event with respect to the trigger particle. In this thesis we will 
concentrate on azimuthal correlations. 

2.1 Analysis Procedure 

We select our trigger particle such that it has transverse momentum, p^, much 
greater than the average pt of the produced particles. These high pt particles may 
predominantly come from jets of hard-scattering partons. The selection of high px 
particles, thus, preferentially triggers on jets. The azimuthal distributions of lower px 
particles, associated particles, in the event are investigated with respect to the trigger 
particles. Wc shall denote the azimuthal angle of the associated particles and trigger 
particles as (j) and ^rrip, respectively. We shall denote the distribution of the lower 
Pt particles with respect to the trigger particle as J2(A0) where A0 — (f) — (pxrig- 

For each trigger particle a distribution is constructed for all particles within a 
given Pt window in the same event. All trigger particles within a given pt range and 
all events within a given centrality are accumulated. The correlation is normalized 
per trigger particle, not per event, because we are interested in quantities on a per- 
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jet basis. The background distribution is constructed through event mixing. 
The associated particles are taken from a different event as the trigger particle to 
construct the A0 distribution. The mixed events must be of the same ccntrality as 
the trigger particle but otherwise no conditions are place on them and will be referred 
to as inclusive events. Our mixed events will take care of dependent effects due to 
detector efficiency to first order. However, in a given event all particles are correlated 
with the reaction plane. This is the flow correlation. The flow correlation between 
the trigger particle and the background particles is lost in event mixing, because there 
is no correlation between the reaction plane of the triggered event and the reaction 
plane of the inclusive event. This correlation is put in by hand using the measured 
elliptic flow {V2) values. In principal, one can mix events with the same (or similar) 
event planes; however, the measured event plane is not equal to the reaction plane 
due to event plane resolution. The constructed mixed event background would need 
to be corrected for the event plane resolution. This is more complicated than the 
current procedure. 

The level of the background created from event mixing is not exact. There are 
two reasons the level is incorrect The first is that choosing events with a trigger 
particle gives us a bias in our event selection towards higher multiplicity events which 
have larger background multiplicity. The background created through event mixing 
is too low to represent the background in the triggered events. This bias is especially 
significant in pp, d+Au, and peripheral Au+Au collisions. Another reason is that 
the total multiplicity used to determine centrality in the triggered events is of the 
underlying event plus the jet. Thus the true background is lower than that from 
mixed events constructed from inclusive events of the same centrality window. This 
effect is usually smaller than the other but can become the dominate effect in central 
Au+Au collisions. To correct for these biases we introduce a normahzation factor 
a. An assumption must be made to the level of the background normalization to 
determine a. The assumption used in our analysis is that the correlation signal is 
zero at = 1 (Zero Yield At 1 or ZYAl) which is the minimum of the correlation 
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signal. This assumption is known as an upper limit on the background since the signal 
is positive definite. A assumption must be made because the true level of background 
is unknown a priori. 



2.1.1 Flow Correction 

In non-central Au+Au collisions, the overlap geometry is not isotropic. This 
anisotropy results in a non-uniform pressure gradient in the initial stage of the col- 
lision, which in turn results in a momentum anisotropy in the final state. This 

anisotropy is characterized by Fourier expansion, 

dN N °° 

— = — [1 + 2 ^ t;„ cosn(0 - ^)] (2.1) 
# 27r ^1 

where (f) and ^ are the azimuthal angles of the particle and the reaction plane, re- 
spectively. At mid-rapidity, the measured directed fiow (fi) is consistent with zero as 
expected due to symmetry. For symmetric collisions, odd orders of Vn are expected 
to be zero at midrapidity. Elliptic flow (^2), however, can be large especially at mid- 
rapidity. For 2-particle correlations the elliptic flow is the only term that signiflcantly 
contributes to our background; all higher orders of Vn are negligible. Keeping terms 
up to order V2 we get the azimuthal distributions of the associated particles and the 
trigger particles respectively: 

^ = ^[l + 2^2cos2(0-*)], (2.2) 

= + 2^''^^ cos 2(0TH. - n- (2.3) 

The correlation functions we are interested in are expressed in terms of A0 = 
(f) — (pTrig- To obtain the flow correction in A0, we integrated over ^, (pTrig, and 4> as, 
dN _ f^^ f^^ f^^ NNrrig I 

[1 + 2^;J"^COs2((/)rrig - *)]5(A(/) - (0 - (t)Trig))d^ d(l)d(l)Trig 

= ^^^"g [1 + 2v2vl"' cos(2A0)] (2.4) 
27r 
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In this thesis, we use the average of the measured f 2 from the modified reaction plane, 
V2{MRP}, and 4-particle cumulant, V2{4}, methods from [IS]. The reaction plane 
method [13] over predicts the V2 due to contributions from non-flow effects (such as 
jets). The modified reaction plane method excludes particles with \Ari\ < 0.5 to 
reduce the non-flow contribution from jets. The 4-particle cumulant method [H] is 
able to suppress the additional contributions from non-flow but under predicts the f 2 
signal in the presence of V2 fluctuations. Figure [2?T| left, shows the V2 values from the 
reaction plane and modified reaction plane methods as a function of Pt for different 
centrality windows. The right panel compares the V2 values for the reaction plane 
and the 4-particle cumulant methods as a function of centrality (along with the values 
from some other methods). Since there are no V2 measurements for the top 5% and 
70-80% for the 4-particle cumulant method we assume f2{4} = V2{MRP}/2 for these 
two centralities. This is a very conservative estimate to account for any systematic 
uncertianties associated with extrapoliation. The ratio V2{4:}/v2{MRP} is plotted in 
figure 12.21 for the measured and estimated values. For our systematic uncertianty on 
the flow correction we will vary the f 2 value used between these two measurements. 
The V2 values for the two measurements are taken from published STAR data [45j. 
The background normalization and flow subtraction proceedue has been previously 
used in Au+Au and pp collisions [25]. 

2.1.2 Other Corrections 

The overall detector efficiency is obtained using a track embedding technique. 
Monte Carlo tracks are placed within real events and then propagated through sim- 
ulations of the detector responses |47j. The probability the embedded track can be 
reconstructed determines the efficiency. The tracks are run through a simulation 
of the detector geometry and then through the tracking software to determine the 
overall efficiency of the detector and the reconstruction algorithm. The efficiency is 
determined as a function of pr using the same track quality and pseudorapidity cuts 



23 




Figure 2.1. Left: Reaction plane and modified reaction plane V2 results 
as a function of pt in different centrality windows. Right: Comparison 
of V2 values using different flow measurement techniques as a function of 
centrality with 'standard' signifying the reaction plane results and 't'2{4}' 
signifying the 4-particle cumulant results. Flow measurements are from 
Au+Au collisions at a/J/vw = 200GeV. Plots are from reference [15] 



as is used in real data. The charged pion efficiencies are used for all tracks because 
the majority of tracks are pions. At low px, the pions dominate; while at high pt, 
the proton contribution is no longer small, but the efficiencies for protons and pions 
are similar [17]. Efficiency corrections are only performed on the associated particle 
tracks only since the correlation functions are normalized per trigger and thus trig- 
ger particle efficiency cancels. Figure 12.31 shows the parameterization to the detector 
efficiency for 0-5% central Au-|-Au collisions. Effiencies for the other centralities of 
Au-|-Au collisions and d-l-Au collsions may be found in the appendix (Fig. IA.II) . The 
fit function is 



Poe pt' 



(2.5) 



24 




Figure 2.2. Ratio of 4-particle cumulant and reaction plane V2- Solid 
points are measured. The line is a fit to a second order polynomial. The 
open points are a conservative estimate of the ratio for the centrality 
bins without 4-particle cumulant V2 measurement. Elliptic flow values are 
from [15] . 



In the analysis the resultant fit function is then evaluated on a particle by particle basis 
and each particle is assigned a weight corresponding to the inverse of the efficiency. 
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Figure 2.3. Detector efficiency for 0-5% most central Au+Au collisions at 
^/sNN = 200 GeV as a function of pr- Errors are statistical. 
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Table 2.1 

Elliptic flow parameterization. For particles below 4 GeV/c V2 = 
-pLPT^e"'^^^/^^^ Above 4 GeV/c the high px f2 is used. The system- 
atic uncertainty on the V2 is given in the last column. 



Centrality 


Pi 


P2 


P3 


P4 


High Pt V2 


Sys. Uncert. 


70-80% 


0.43383 


1.587 


0.514 


0.568 


0.19265 


±33% 


60-70% 


0.26470 


1.410 


1.376 


0.809 


0.23667 


±12% 


50-60% 


0.19412 


1.258 


2.241 


0.981 


0.19747 


±9% 


40-50% 


0.20927 


1.303 


2.002 


0.981 


0.19267 


±7% 


30-40% 


0.16241 


1.208 


2.750 


1.098 


0.18034 


±6% 


20-30% 


0.14261 


1.213 


2.846 


1.027 


0.17471 


±6% 


10-20% 


0.12437 


1.308 


2.326 


0.883 


0.13340 


±6% 


5-10% 


0.08216 


1.341 


2.357 


0.924 


0.09434 


±12% 


0-5% 


0.03727 


1.273 


3.133 


1.352 


0.04797 


±33% 



2.2 d+Au Collisions 

Two-particle jet-like correlations have been previously studied in Au±Au and pp 
collisions [25]. These results have shown a broadened and enhanced away-side peak 
at low associated px and a suppressed away-side peak at high associated pt in central 
Au±Au collisions with respect to pp collisions. Two-particle azimuthal correlations 
in d±Au collisions at high px have been studied before [23]. In this thesis, we analyze 
d±Au collisions with associated particles going down to low pt- These collisions are 
interesting because they will allow us to disentangle initial state nuclear effects from 
other effects seen in Au±Au collisions when compared to pp collisions. 
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2.2.1 Results and Discussions 

We study these correlations as a function of A?] and A0, the A0 correlations are 
studied in greater detail. Figure 12.41 shows the raw correlation functions in Arj on 
the left and A(f) on the right in the top panels with the open red signals. The Arj 
correlations are shown only for particles on the near side (particles with |A0| < 1). 
The black lines in the top panels show the backgrounds from event mixing. Since this 
is d+Au no f 2 correlation is present so our A0 background is just a flat line. The effect 
of the TPC sector boundaries is too small to be visible on this plot. The triangle-like 
shape to the At] background is a result of the detector 2-particle acceptance. The 
triangle-like shape is not symmetric about At] = because d+Au is an asymmetric 
system. The trigger particle is constrained in \Ai]\ < 0.7 while the associated particle 
are in |A?7| < 1.0. The associated particles are given the larger At] range so that most 
of the near-side jet-cone will be within our acceptance. The bottom left and right 
panels show the background subtracted signal for At] correlations and A0 correlations 
in red. The background we subtracted does not quite match the magnitude of the true 
background, as seen in the upper right and lower left panels. There are 2 competing 
reasons for this mismatch. The first (and the dominate one in d+Au collisions) is 
when we trigger on a high pt particle we bias our events towards higher multiplicity. 
This leads to the background being too low. The second arises from the centrality 
selection by cutting on the total reference multiplicity. Since the total reference 
multiplicity includes jet-like correlated particles, the true background level is smaller 
than the mixed-event multiplicity. We scale our mixed event background such that 
after subtraction the jet-like correlation signal is zero at A0 = 1. This is done over 
a fixed range |A0 ± 1| < 0.2. The same scaling factor has been used for A0 and 
At] correlations. The final correlation functions are then obtained by subtracting the 
normalized background. The blue points in the Arj correlation have this additional 
mixed event background subtracted. 
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Figure 2.4. Analysis plots for minimum bias d+Au at a/s/vjv = 200 GeV/ c. 
Top Left: raw signal in At] (open red points) and background from mixed 
events (black line). Bottom Left: background subtracted signal in Ai] (red 
points), background subtracted signal with additional normalized back- 
ground subtracted (blue points), and residual background obtained from 
ZYAl (black line). Top Right: raw signal in A0 (open red points) and 
background from mixed events (black line). Bottom Right: background 
subtracted signal in A0 (red points) and residual background obtained 
from ZYAl. The pt ranges of the trigger particles is 4 < < 6 GeV/c 
and that of the associated particles is 0.15 < < 4.0 GeV/c. Errors are 
statistical. 
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The A0 signal with full background subtraction for d+Au collisions is shown if 
Figure 12.51 in the red symbols for 4 < p^"^ < 6 GeV/c and 0.15 < p^"'^"^ < 4 
GeV/c. The left panel has been fit to Gaussians with centroids fixed at and tt 
for comparison. The right panel shows a comparison to the published pp correlation 
(shown in blue points). The d+Au data has been rebinned in this plot. The d+Au 
and pp correlations are consistent which indicates that there is not significant initial 
nuclear effect in d+Au collisions. It implies that the modification of the away-side 
jet-like correlation shape observed in central Au+Au collisions relative to pp collisions 
is not an initial nuclear effect. 
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Figure 2.5. Background subtracted distributions for d+Au collisions 
at v^s^ = 200 GeV/c (red) for 4 < p^"^ < 6 and 0.15 < < 4 

GeV/c. Left: Fit with two Gaussians with centroids fixed at and tt. 
Right: Compared with pp [2^ (blue). Errors are statistical. 



We study the correlation functions versus associated particle pt and extract the 
correlated yields versus associated px- Figure shows the associated particle spec- 
tra. The left column of plots show the near-side spectra and the away-side spectra 
on the top and bottom, respectively. The spectra are shown for trigger particles of 
3 < < 4 GeV/c, 4 < < 6 GeV/c, and 6 < pr < 10 GeV/c. The spectra 
become harder as the trigger pt increases on both the near-side and away-side. The 
right column of plots shows the Z-spectra where Z = p^'^'^°^ / p^"^ . The spectra are 
shown this way to present it in a format that is closer to the fragmentation function 
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which is in terms of z = p^s'^o'^ jpJ^^ ^ where p"^*^* is the parton momentum (or energy) 
and and p^**""^ is the fragment momentum projection along the parton momentum 
direction. The Z spectra are harder for the lower px triggers, as seen in the figure. 
The real fragmentation function in z is independent of the jet energy. The reason for 
the change as a function of trigger pj- is that the trigger particle momentum becomes 
a better proxy for the parton momentum as we increase Pt-, our triggering on a high 
Pt particle biases us towards jets where more of the parton energy is contained in the 
trigger momentum this bias increase with trigger px- Figure 12.71 shows the near-side 
and away-side spectra for both d-l-Au and pp collisions. The spectra are constant 
between the two collision systems showing no significant initial nuclear effects. 



Near Side Spectra | I Mear Side Spectra | 




Figure 2.6. Associated particle spectra in d-l-Au collisions at ^/snn = 200 
GeV/c. Left: spectra. Right: Z = spectra. Top: Near-side 

spectra. Bottom: Away-side spectra. Errors are statistical. 
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Figure 2.7. Comparison of d+Au and pp spectra at ^/sJ^ = 200 GeV/c 
for 4 < py"^ < 6 GeV/c. Both the near-side and away-side spectra are 
shown. The pp spectra is from [25] . Errors are statistical. 



The width of the jet-hke correlations is also interesting. Figure [22] (left) shows the 
RMS width as a function of associated particle pt for both pp and d+Au collisions 
for 4 < p^"^ < 6 GeV/c. Widths are shown in both r] and for the near-side and for 
on the away-side. The width of the rj distribution on the away-side is not shown 
because the away-side is flat within our detector acceptance. This is because the 
near-away and away-sides partons are not correlated in rj due to the colliding partons 
carrying different fractions of the longitudinal momentum to the collision giving the 
jet an overall longitudinal momenta of the colliding nucleons. The d-|-Au widths show 
no significant changes from the width in pp. The right panel shows the width as a 
function of trigger particle pt for associated particles of 0.15 < px < 3.0 GeV/c. The 
open symbols are the widths from Gaussian fits. The filled symbols are from RMS. 
The widths are consistent for the two measurement techniques because the peaks are 
very Gaussian. In both figures, the widths are similar in rj and on the near-side. 
This shows that the near-side emission is within a cone. The near-side should not 
have same width as what would be calculated in fragmentation for widths respect 
to the jet-axis. It should be broadened due to the use of the trigger particle as a 
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proxy for the jet-axis. The away-side has additional broadening due to two coupled 
reasons. One is (initial parton transverse momentum). This initial transverse 
momentum results in the near-side and away-side not being quite back-to-back in (p 
for a given di-jet. Upon averaging over many di-jets this results in a broadening on 
the away-side. The other reason is trigger bias, where the trigger particle tends to 
select the jet with higher transverse momentum, also due to the fact that the partons 
carry some initial transverse momentum. The backside partner jet has lower energy 
and thus the fragmentation cone is wider. In principle, one should be able to extract 
valuable information about the initial kx broadening from the width measurements 
in pp and d+Au within a model framework. One interesting note is that one would 
expect the fcy broadening to be larger in d+Au than in pp because of initial multiple 
scatting in d+Au, however, our pp and d+Au comparison result does not seem to 
show this is a significant effect. 



dAu MinBias 0<pt assoc<3 | 




Figure 2.8. Jet-like correlation peak widths. Left: RMS widths as a 
function of associated pt- Open symbols are for pp [25] and filled symbols 
are for d+Au ^/sNN = 200 GeV/c. Circles are for the away-side width in 
0. Squares are for the near-side width in 0. Triangles are for the near-side 
width in rj. Right: RMS width (filled) and Gaussian fit sigma (open) as 
a function of trigger px- Errors are statistical. 
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2.2.2 Summary 

Jet-like correlations have been studied in d+Au collisions at ^/snn = 200 GeV/c. 
The correlations are consistent with jet-like correlations in pp collisions. This implies 
that the modification to the away-side correlation observed in central Au+Au col- 
lisions is not an initial nuclear effect. This also justifies the use of d+Au data for 
comparison to Au+Au collisions instead of pp collisions. The benefit gained from 
this is the greatly increased statistics for d+Au collisions compared to pp. Associated 
particle spectra and correlation widths are extracted for d+Au collisions and can be 
used for comparison to theoretical calculations. 

2.3 Au+Au Collisions 

Two-particle correlations have already been studied in Au+Au collisions at ^/snn = 
200 GeV at STAR [IHIES]. The publications used Au+Au data from the second year 
of RHIC running. This thesis extends the measurements using data from the fourth 
year of RHIC running which has about an order of magnitude increased statistics. 
In this thesis, we used the same trigger and associated pt cuts as are used in the 
PHENIX publication ^8] to check for consistency between the two experiments. 

The background subtraction method in Au+Au collisions is the same as that for 
the d+Au collisions except for the fiow subtraction. The elliptic fiow has been added 
to the mixed events pairwise by hand as discussed in section 3.1.1. The top panel 
of Figure 12.91 shows the raw 2-particle correlation function in in red. In black 
are the mixed events, where the associated particles are taken from different events 
of the same centrality as the trigger particle. In blue are the mixed events with the 
elliptic fiow correction added in pairwise. The bumps and dips apparent in all three 
distributions are due to the TPC sector boundaries. The bottom panel shows the 
background subtracted correlation function. This is the raw 2-particle correlation 
function minus the normalized mixed event background with elliptic fiow. The mixed 
event with elliptic fiow background is normalized such that the background subtracted 
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Figure 2.9. Two-particle correlation functions for 2.5 < p-p^*^ < 4.0 GeV/c 
and 1 < p^''"^ < 2.5 GeV/c in 10-20% Au+Au collisions at ^/s^ = 200 
GeV/c. Top: Raw correlation function (red). Background from event 
mixing (black). Background from event mixing with elliptic flow added 
in pairwise by hand (blue). Bottom: Background subtraced 2-particle 
correlation function. The background is normalized by ZYAl before sub- 
traction. Errors are statistical. Plots for all of the Au+Au centrality bins 
can be found in the appendix. 



correlation function is zero in the region |A0 ± 1| < 0.2 [25] as was done in d+Au 
collisions. 

2.3.1 Comparison Between STAR and PHENIX 

Figure 12.101 shows the background subtracted 2-particle correlation functions in 
Au+Au collisions at ^/sNN = 200 GeV/c using the high statistics year 4 data. The 
correlation functions are for 2.5 < Pt"'' < 4.0 GeV/c and 1.0 < p^''°^ < 2.5 GeV/c. 
The left set of panels shows our results and the right panels show the published 



34 



STAR Preliminary 
0-5% 



STAR Preliminary 

5-10% 




A(|) (radians) 




0.4 
0.35 

0.3 
25 

02 
0.15 

0.1 
05 



0.35 
03 

25 
02 

0.15 
0.1 

05 






0.5 1 1.5 2 2.5 3 
A i|i (rad) 



; 1.5 2 2.5 3 
A([i (rad) 



Figure 2.10. Background subtracted correlation functions in Au+Au 
collisions at = 200 GeV/c for 2.5 < p^"^ < 4.0 GeV/c and 

1.0 < < 2.5 GeV/c. Left: Our results with Iri"^^^"^] < 1.0. Right: 

PHENIX results from [48j with < 0.35. Error bars are statistical 

errors. Histograms represent 1 sigma systematic uncertainty on the ellip- 
tic flow, except that the top most histogram on the PHENIX results is -2 
sigma systematic uncertainty. Shaped bands on PHENIX results are the 
normalization uncertainty. 



PHENIX results [IH]. The magnitude of the peaks is higher in our results than 
the PHENIX results due to our larger acceptance of associated particles in r]; our 
results have |r/^^«°c| < 1.0 while the PHENIX results are for < 0.35. This 

increases the magnitude of our away-side peaks by 1.0/0.35 since the away-side is 
evenly distributed in rj [25]. Our near-side peak is also larger because of two reasons. 
One is that the near-side peak is broad and our larger rj acceptance catches more 
of the associated particles. The other is that we include more of the long range t] 
correlations (the ridge) that have been observed in 2-particle A0 — Arj correlation 
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Figure 2.11. Away-side correlation amplitudes in the hump (7r/3 < |A</) — 
7r| < 47r/9) and center (|A0 — vr| < vr/9) as a function of centrality, 
{Npartf^y^^ shown in blue and red respectively. Our data is shown in 
filled symbols and has been scaled down by a factor of 0.35 to account 
for the rj acceptance. The PHENIX data is shown in open symbols. The 
bars shown the statistical errors and the systematic errors due to flow 
are shown in shaded bands and caps for our data and the PHENIX data, 
respectively. 



We are most interested in comparing the shape of the away-side correlation be- 
tween the two experiments. Our results show a dip at = tt in the three most 
central bins; however, this result is not systematically significant given the system- 
atic errors (shown by the black histograms) due to the uncertainty on the elliptic 
flow result. In our results, the dip is the strongest in the 0-5% most central bin as 
would be expected if the dip is a medium-induced effect. In the PHENIX results, 
the dip at A0 = vr is stronger and is systematicly significant with their uncertainty 
on elliptic flow (three histograms are ±lcr and —2a). The dip at tt in the PHENIX 
results is the strongest in mid-central collisions where the elliptic flow uncertianities 
are the largest. This is illistrated by Fig. l2.10l where the centrality dependences of the 
hump (7r/3 < \A(f) — 7r| < 47r/9) and center (| A0 — 7r| < tt/Q) are shown for both our 
results and the PHENIX results. The trends of the centrality dependences suggests 
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that the discrepancy between our results and the PHENIX results may come from the 
systematics in the flow subtraction. This figure also shows that the signal strength in 
the central region does not drop as centrality increases. The double peaked structure 
in central Au+Au collisions results from an increased yield in the hump region. 
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Figure 2.12. Background subtracted correlation functions in Au+Au 
collisions at = 200 GeV/c for 2.5 < p^"^ < 4.0 GeV/c and 

1.0 < p^''°^ < 2.5 GeV/c. Left: Our results with |r^^^«°'=| < 1.0. Right: 
Our results with |7^^*«°'=| < 0.35. Error bars are statistical. Histograms 
represent 1 sigma systematic uncertainty on the elhptic flow. 



A cross check has been done to see if there is any effect in the away-side shape due 
to the r] acceptance. Figure [2.121 shows a comparison of our results with < i.Q 

and < 0.35. No significant difference on the away-side shape is seen. 
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Figure 2.13. (a) Proton to tt"*" and (b) anti-proton to vr" ratios for Au+Au 
collisions in two centralities and d+Au collisions at y^s^N = 200 GeV/c 
from ^19j . 



2.4 Identified Trigger Particle Correlations 

An enhanced baryon to meson ratio at intermediate pt has been observed in 
central Au+Au collisions with respect to peripheral Au+Au, d+Au and pp colli- 
sions [5n[52] . Figure [2.131 shows the proton(anti-proton) to 7r+ (vr^) in the left (right) 
panel for central Au+Au, peripheral Au+Au and d+Au. This effect can be explained 
by coalescence and recombination models [531[Sll|35j. In this model, the momentum 
comes from the sum of the quark momenta used to make the particle. Since baryons 
are composed of three quarks and mesons of two the baryons are shifted to higher 
momentum than the mesons resulting in an enhanced baryon over meson ratio at 
intermediate pt- However, simple coalescence and recombination models do not con- 
tain any angular correlations, thus one would expect lower per trigger correlation 
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strength for intermediate pt triggered baryons then triggered mesons. In this section, 
we identify intermediate pt baryons and mesons and using them for trigger particles 
to do jet-correlation studies. We attempt to identify and study any differences in the 
associated particle distributions to hopefully shed further light on the baryon/meson 
puzzle. 



2.4.1 Particle Identification in TPC at High Pt 




Figure 2.14. Ionization energy loss plotted as /o5fio(dE/dx) vs. logioip). 
The lines represent ±0" (±8% in relative dE/dx) bands for protons, kaon, 
electrons, and pions. 170 stands for Bichsel's prediction for 30% truncated 
mean for dE/dx. Figure is from [56] . 



The STAR time projection chamber (TPC) was designed for identification of pro- 
tons, anti-protons and charged pion and kaons at low pt using energy loss (dE/dx). 
At around 1.5-2.0 GeV/c the energy loss is similar enough for all three particles that 
they cannot be distinguished. However, particle identification can be performed for 
higher px particles. This is due to the relativistic rise of dE/dx. The dE/dx distribu- 
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tions for different particle types is shown in figure [2?T^ as a function of the magnitude 
of momentum. Pions can be distinguished from kaons and protons at around 2.5 
GeV/c and above. Protons and kaons can be distinguished from each other starting 
at about 3.0 GeV/c. 
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Figure 2.15. A^cr^ distribution for 3 < < 3.5 GeV/c and |?7| < 0.35 
from 0-5% Au+Au collsions at ^/sNN = 200 GeV/c. Positive particles are 
plotted as NaT, plus 10, negative particles are plotted as -Na-,, minus 10. 
Curves are from fit, red is pion, green is kaon, blue is proton, and black 
is the sum of the three. Errors are statistical. Plots for other bins can be 
found in the appendix. 



To separate the particles Na-,, is used, where Na-,, is the (dE/dx-predicted n 
dE/dx) /resolution. With this variable, if the data is well calibrated, the pions will 
form a Gaussian peak centered at zero with a a width of 1. Figure 12.151 shows an 
example Na.,, distribution. In this figure, positive and negative particles are separated 
by plotting ±A^(Tvr ± 10 for positive and negative particles, respectively. It is beneficial 
to bin the dE/dx distribution in rj in addition to pt and centrality. This is because 
the energy loss is a function of the magnitude of the total velocity and with too large 
of an rj bin we are sampling a large range of total momentum magnitude within one 
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Pt bin, resulting in a reduced resolution. We can than fit the Na.^ distributions 
to six Gaussians, one each for 7r+, vr^, K~^, K~~ , p and p. However, this fit needs 
additional constraints due to overlap of the Gaussians, as seen in Fig. I2.15[ One 
set of constraints is the requirement of the same width and centroid (adjusted for 
the way it is plotted) for the positive and negative particles. Since the protons and 
kaons are not well separated in this region from just energy loss more constraints are 
required to extract the protons. We can use the Bethe-Block formula to calculate 
the energy loss. This can be used to determine the kaon centroid position relative 
to the pion and proton. From this calculation the kaon peak position is determined 
to be approximately 1/3 of the distance between the proton and pion peaks. The 
calculated centroids cannot be simply fixed in the fit because the dE/dx calibration 
is not perfect. It may also be necessary to use the measured yield for the charged 
kaon yield, with the assumption that has the same yield as and K~ . An 
example fit is shown in Fig. 12.151 This method of fitting the dE/dx distributions has 
been previously used to extract high pt identified particle spectra [T9] . 

Cuts can then be placed on the NaT^ distribution to extract enriched samples of 
protons, anti-protons, and charged pions. We can use our fit results to determine the 
purity of the selected particles from the cuts. If there is large overlap of the Gaussians, 
we may have to sacrifice purity for statistics. Figure 12.161 shows an example Na^^ 
distribution for the positive and negative particles with lines designating where cuts 
are to obtain 50% purity for p/p and 95% purity for tt^/tt^. Particles with Nct, 
below the left line are taken to be protons and particles above the red line are taken 
to be vr+s with the given purities. The purity is determined by integrating the curve, 
within the particle identification cut, obtained for the fit for the particlar particle 
type and dividing by the intergral of the curves for all particle types. 

The position of these cuts can be improved upon then by fitting the cuts to 
a smooth curve as a function of px for each rj bin. An example fit is shown in 
Figures I2.17[ The cuts were computed out to higher px but due to statistics only 
trigger particle of 3 < < 4 GeV/c were used. The fits are to second degree 
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Figure 2.16. Na.,, distribution for 3 < pt < 3.5 GeV/c and \rj\ < 0.35 from 
0-5% Au+Au coUsions at a/s/vat = 200 GeV/c. Left: Positive particles. 
Right: Negative particles. Curves are from fit, red is pion, green is kaon, 
blue is proton, and black is the sum of the other three. The vertical lines 
show the position of the cuts. Na.,, greater than red line gives 95% purity 
of charged pions. Na.„ less than blue hne give 50% purity of proton/anti- 
protons. Errors are statistical. Plots for other bins can be found in the 
appendix. 
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Figure 2.17. Fits to the cuts on A^a^ as a function of pt for vr^ with 95% 
purity. Plot is for Au+Au 0-5% central collisions at ^Jsnn = 200 GeV/c 
with |?7| < 0.35. Errors are statistical. 
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polynomials. Positive particles with Na.,^ greater than the vr"*" fit and negative particles 
with A^cTtt greater than the vr" fit are considered charged pions with 95% purity. 
Positive particles with Na^^ less then the p fit and negative particle with NaT^ less 
than the p fit are considered protons/anti-protons with 50% purity. 

2.4.2 Results and Discussions 

Figure 12.181 shows identified trigger particle correlations in Au+Au collisions for 
3 < py"^ < 4 GeV/c and 0.15 < p^^'^°'^ < 3 GeV/c. In this region, the correlations are 
dominated by the lower associated particles. The raw signals are shown in the top 
panel of each set of three panels. The difference in the background levels seen for the 
raw signal is due to different trigger biases for different types of triggers. The second 
panel shows the raw signal minus mixed events. The mixed events are normalized 
such that the signal is ZYAl. Once the difference level of the background is removed 
the correlations look similar. Part of the correlation is not jet-like but is due to 
the correlation with the reaction plane. The bottom panel shows the correlations 
with the elliptic flow modulated mixed events removed. The flow is added into the 
mixed events pair-wise. The flow is not the charged particle flow as used above, 
but is scaled such that the baryon/meson flow is 3/2. This contains the assumption 
of quark number scaling to the elliptic flow which has been seen in data up to 2 
GeV/c [57]. The yields obtained from the Gaussian fits in the particle identification 
procedure are used to constrain the relative number of baryons and mesons. After 
elliptic flow subtraction the jet-like correlations show no significant dependence on 
trigger particle type for all centrality bins. Figure 12.191 shows the identified trigger 
particle correlations in Au+Au collisions for 3 < p^"^ < 4 and 2 < p^^^°^ < 3 GeV/c. 
These higher px associated results are similar to the lower pt associated results in that 
background subtracted jet-like correlations show no significant dependence on trigger 
particle type. Figure 12.201 shows the identified trigger particle correlations in d+Au 
collisions for 3 < p^"^ < 4 and 0.15 < p^''°^ < 4 (left) and 2 < p^''°^ < 4 GeV/c 
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Figure 2.18. Identified trigger particle correlations in Au+Au collisions 
at -y/s/vTV = 200 GeV/c. Each group of 3 plots are (from top to bottom) 
raw signals, raw signal minus mixed events, and raw signal minus mixed 
events and V2- Groups from left to right, top to bottom are centralities 50- 
80%, 30-50%, 20-30%, 10-20%, 5-10%, and 0-5%. All plots are for Au+Au 
collisions with 3 < p?"^ < 4 GeV/c and 0.15 < p^''°^ < 3 GeV/c. Errors 
are statistical. 
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(right). For the d+Au plots, the p and p and the vr"*" and vr" have been combined 
due to low statistics. Again for the identified trigger particle correlations in d+Au 
no significant differences are seen for the triggered baryons and mesons. 
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Figure 2.20. Identified trigger particle correlations in d+Au collisions at 



= 200 GeV/c. All plots are for 3 < p^^ < 4. Left: 0.15 < p^'""' < 
3. Right: 2 < p^'^^°^ < 3. Top: Raw signals. Bottom: Raw signals minus 
normalized mixed events. Errors are statistical. 



Fi gures [2 . 2 1 1 and 12 . 2 2 1 show the near-side and away-side associated particle spectra 
for triggered protons/anti-protons, charged pions, and all charged particles. The 
associated particle spectra for both triggered protons/anti-protons and charged pions 
are consistent with the associated spectra for unidentified triggers. This is true for 
all centralities in Au+Au collisions. 

2.4.3 Summary 

Jet-like correlations have been studied with identified trigger particles at ^/sNN = 
200 GeV/c. This study is motivated by the observed baryon/meson puzzle (the large 
baryon/meson ratio at intermediate px)- If more baryons are indeed formed from 
coalescence of thermal constituent quarks, then the per-trigger normalized jet-like 



46 




" 0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 




Figure 2.21. Near side spectra for identified protons/anti-protons (blue), 
charged pions (red) and all charged particles (black) in Au+Au collisions 
at y/SNN = 200 GeV/c. Panels are d+Au minimum bias and Au+Au 
centrality bins 50-80%, 30-50%, 20-30%, 10-20%, 5-10%, and 0-5% from 
left to right and top to bottom. Errors are statistical. 



correlation will be reduced for trigger baryons relative to that for trigger mesons. 
This analysis is done for intermediate pt triggers for two reasons. First, we are 
statistics limited to explore the p^^^^ space in more detail. Second, this is the region 
of interest where the baryon/meson ratio peaks. No significant difference is seen 
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Figure 2.22. Away side spectra for identified protons/anti-protons (blue), 
charged pions (red) and all charged particles (black) in Au+Au collisions 
at y/SNN = 200 GeV/c. Panels are d+Au minimum bias and Au+Au 
centrality bins 50-80%, 30-50%, 20-30%, 10-20%, 5-10%, and 0-5% from 
left to right and top to bottom. Errors are statistical. 



between jet-like correlations with identified baryons {p and p) and mesons (vr" and 
TT^). The increased baryon/meson ratio does not seem to have a significant effect on 
the associated particle production. 
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3. THREE-PARTICLE JET-LIKE CORRELATIONS 

3.1 Introduction and Motivation 

Two-particle azimuthal correlations that are reported in Chapter 3 of this thesis 
and elsewhere [251H8] have shown a modified away-side structure in central Au-|-Au 
collisions, with respect to pp, d-l-Au, and peripheral Au+Au. The away-side distri- 
bution is broadened or even double-humped in central Au+Au collisions. Different 
physics mechanisms have been suggested to explain this modification including: large 
angle gluon radiation [58|[59]. jets defiected by radial fiow [60] or preferential selection 
of particles due to path- length dependent energy loss [^I], hydrodynamic conical fiow 
generated by Mach-cone shock waves [62l[63l[6ll[65l[66l[67l[68l[69], and Cerenkov gluon 
radiation [701171] . Three-particle correlations can be used to differentiate the mech- 
anisms with conical emission, Mach-cone and Cerenkov gluon radiation, from other 
physics mechanisms. 

Figure 13.11 is a cartoon view of expected 3-particle correlation signals for different 
scenarios. The top panel shows unmodified back-to-back di-jets. In A^i — A02 space, 
, where A0j = (pi — 4>Trig, the 3-particle correlation will have four peaks, one at (0, 0) 
for instances where both associated particles are on the near-side, one at (vr, vr) for 
instances where both associated particles are on the away-side and peaks at (0, vr) 
and (vr, 0) for instances where one of the associated particles is on the near-side and 
the other is on the away-side. The middle panel shows what happens in situations 
where the away-side jet is modified such that particles only come out to one side of 
vr. In a single event, this would result in an away-side peak on-diagonal but displaced 
from vr; however we average over many events resulting in an on-diagonal structure. 
The structure stays on-diagonal because even though the away-side particles are not 
at vr they are still close together. This scenario can come about due to physics 
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mechanisms such as: jets deflected by radial flow or path-length dependent energy 
loss. The bottom panel shows the expected cartoon results for conical emission. For 
conical emission particles can be off-set to both sides of tt. When both associated 
particle are off-set to the same side of vr we get on-diagonal peaks again. If we take 
particles that are off-set to opposite sides of vr the particles are no long close together 
and we get peaks on the off-diagonal. This gives us four peaks on the away-side. The 
two off-diagonal peaks are only expected for conical emission and therefore will be 
our signature for conical emission. 

3.2 Analysis Procedure 

Results from this analysis have been publicly shown by STAR [721I731I7H [751176] . 
This chapter presents the details the analysis. The details of the analysis techniques 
have been described in [77]. The analysis assumes that an event is composed, besides 
the trigger particle, of two components, particles that are correlated with the trigger 
particle (other than the correlation due to anisotropic flow) and background particles 
from the bulk medium (only correlated with the trigger particle due to anisotropic 
flow). Suppose the number of jet-correlated particles is Nj^t and the number of 
background particles is Nbkgd- The total number of particles is then N = Njet + Nbkgd- 
The number of particle pairs is N{N—1) and is composed of three parts: the number of 
background pairs Nbkgd{Nbkgd — 1), the number of jet-correlated pairs Njet{Njet — 1), 
and the number of cross pairs NjetNbkgd + N^kgdNjet = '^NjetNbkgd- Since we are 
interested in the angular correlation of the jet-correlated pairs which is an unknown, 
we must subtract the angular correlations of the background pairs and the cross pairs 
that are known. The angular correlations of the background pairs and the cross 
pairs can be obtained because they are uncorrelated (except for the anisotropic flow 
correlation). Therefore, by subtraction from the raw 3-particle correlation the angular 
correlations of the background pairs and the cross pairs, one can obtain the genuine 
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Figure 3.1. Cartoon representations of the expected results for 3-particle 
azimuthal correlations for different physics mechanisms. Top: Unmodified 
back-to-back di-jets. Middle: Jets deflected by radial flow or path length 
dependent energy loss. Bottom: Conical emission. 



3-particle jet-correlation. Lets introduce some equations to restate this more formally. 
For 2-particle correlations: 



J2(A0) = dNld/\(p, 



51 



S2(A(/)) = dNbkgd/dA(j), 
J2(A0) = dNjet/dA<i), 

J2(A0) = J2(A0)-52(A0), (3.1) 

where A0 is the associated particle azimuthal angle relative to that of the trigger 
particle. J2 is a notation for the raw azimuthal distribution relative to the trigger 
particle azimuth for the entire event. B2 is a notation for the azimuthal distribution 
of the underlying background particles relative to the trigger particle azimuth. J2 is 
the azimuthal distribution of the associated jet particles relative to the trigger particle 
azimuth. 

The raw 3-particle azimuthal correlation between the trigger particle and two 
other particles in this notation is: 

J3(A(/)i,A02) = J2(A0i) ® J2(A02) 

= [J2(A0i) + 52(A0i)]® [J2(A02) + 52(A02)] (3.2) 

where A<^B is not a simple product of A and B due to correlations, A0i — (f)i — (f)trig, 
and A02 = 02 — <Ptrig- The raw 3-particle correlation can be divided into three parts: 
the 3-particle jet- correlations we are interested in, 

J3(A0i, A02) = J2(A(/)i) (8) J2(A02), (3.3) 

and two background terms. One of the combinatoric backgrounds is between an 
associated jet particle and a background particle, 

J2®52 = J2(A0i) ®52(A02) +52(A0i) ® J2(A02) 

= J2(A0i)S2(A02)+52(A0i)J2(A02), (3.4) 

which we term the hard-soft background. Since the jet-correlated particle and the 
background particle are uncorrelated, the hard-soft background can be obtained from 
the simple product of the 2-particle jet-correlation function J2 and the background B2. 
The other background is the combinatoric background between the two background 
particles, 

53(A0i, A02) = 52(A0i) ® 52(A02), (3.5) 



52 



which we term the soft-soft background. The soft-soft background can be obtained 
from inclusive events and will be discussed in greater detail below. The final 3-particle 
correlation can thus be determined by, 

J3(A0i, A02) = J3(A(/)i, A02)-[J2(A0i)S2(A02)+S2(A0i)J2(A(/)2)]-S3(A0^ 

(3.6) 



3.2.1 Flow Correction 

In heavy-ion collisions, particle emission is correlated to the reaction plane due to 
hydrodynamic collective flow of the bulk medium and the anisotropic overlap region 
between the coUiding nuclei. We have already examined this flow correlation in 2- 
particle correlation studies up to the second order harmonic. Because the signal to 
noise ratio in 3-particle correlations is small, the fourth order harmonics of the flow 
correlation cannot be neglected. The flow correlation, expressed in harmonics up to 
fourth order, is given by, 
dN N 

— = — [1 + 2^2 cos 2(0 - + 2v4 cos 4(0 - *)] (3.7) 

d(p 2% 

where (p is the azimuthal angle of the particle and ^ is the azimuthal angle of the 
reaction plane. Due to symmetry the flrst and third (and all other odd) harmonic 
terms vanish at mid-rapidity for symmetric collisions systems. The distribution of 
the trigger-background particle pairs is given by 

_ ^[i + 2^*"^cos2(0,,,,-*) + 2^;r^cos4(0,,,,-*)] 



d(j)trigd(j)d{/\4>) 27r 

-[1 + 2^2 cos 2(0 - 2^) + 2vi cos(40 - 4*)] 



Nbkgd f 

•^(5(A0-(0-0iHs)) (3.8) 

for a given ^, (j)trig-i and 0. Nt^ig is the number of trigger particles and N^kg^ is the 
total number of background particles. Integrating over ^ , (ptrig, and we obtain, 
dN _ /•2- f^- d^N 
d(A0) - io io io ^^^'^^'^*"^d*d0d0,H. 

= ^'"sNbkgd ^^ ^ ^^trig^^ cos(2A0) + 2vt%4 cos(4A0)]. (3.9) 
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When normalized per trigger particle, the anisotropic flow modulation is then, 

52(A0) = Bi[l + 2v2'^V2 cos(2A0) + 2v'r^V4 cos(4A0)] (3.10) 

where Bi = N^kgd/ 2tt is the average background density level. 

The distribution of the trigger-background-background triplets is given by. 



■[1 + 2^^^ cos(20i - 2^) + 2vl cos(40i - 4^)] 



....^...^..^..(A^OcilA^.) - ^[l+2.--s(2,,„,-2^) + 2.-cos(4^.„,-4v.)] 

^1 r. , . fl) 



27r' 

■ — [1 + 2vi^^ cos(202 - 2^) + 2vl cos(402 - 4^)] 
27r 

■ — (5(A0i - (01 - ^tr^g))S{A^2 " (02 " ^trig)) (3.11) 

for a given set of 4>trig, (pi, and 02. After integrating over 4>trig, 0i, and 02 we 
obtain, 

d^N /-s^ r^^ /-s^ /•2^ , , . , d^N 



d^!d(j)-, 



rf(A0i)rf(A02) 7o 7o Jo Jo '^'^^rf^rf0jri<;ci0id02rf(A0i)rf(A02) 
= ^^^^^|^[l + 2.^^^^^^^^ 

+24"^4^^ cos(2A02) + 24^^4^^ cos(2A0i - 2A02) 

+2T;f^4'^ cos(4A0i) + 2vi"'vP cos(4A02) 

+2t;l^^t;f ^ cos(4A0i - 4A02) + 2v2'%^\f^ cos(2A0i - 4A02) 

+2i;*"^t;f ^4'^ cos(4A0i - 2A02) + 2^')^;^ ^ cos(2A0i + 2A02)]. 

(3.12) 

Here f ^ is the n*^ order harmonic of the particle i (where i can be the trigger, or either 
one of the two associated particles). In Eqn. 13.121 we have taken the two background 
particles from two separate sets with multiplicities A'^i and N2 giving us N1N2 pairs. 
For particles from the same set that becomes Nbkgd{Nbkgd — ^)- Normalized per trigger 
particle, the anisotropic flow modulation is then given by, 

S3(A</'i,A</'2) = ^''^'^^'''^ + 2vr'vi'^ cos(2A0O 
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^-2v';'Hf cos(2A(/.2) + Iv^i^vf cos(2A0i - 2A02) 
^-2vT'vf cos(4A0i) + 2vT'vf cos(4A02) 
+24^)ff ^ cos(4A0i - 4A02) + 2t;2"^4^^ff ^ cos(2A0i - 4A02) 
^2v';'%fvi^ cos(4A0i - 2A02) + Iv^^^vi^'r cos(2A0i + 2A02)] 

(3.13) 

As in 2-particle correlations we use the average from the reaction plane and 
4-particle cumulant measurements. Since we are now using the fourth order har- 
monics we also need f4 values. Figure [22] shows the measured v^lv\ [45]. We fit 
the measurement of v^lv\ to a constant in the range of IjPTjS GeV/c, yielding 
^4 = \.\^v\. 
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Figure 3.2. Results for f4/f2 [45j. A fit to a constant in the range of 
1 < < 2 GeV/c is used to parameterize f4 in terms of v\. The black 
line shows a fit to a constant over the entire "Pt range. The red curves are 
from blast wave fits. 
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3.2.2 Acceptance and Efficiency Corrections 

All associated particles are corrected for the overall 0-averaged detector efficiency 
as is discussed in 3.1.2. The dependent efficiency (acceptance) is accounted for 
separately. The (f) dependence of the efficiency can be corrected for to first order by 
mixed events. This is sufficient in 2-particle correlations (i.e. the 2-particle acceptance 
correction) but not for 3-particle correlation due to the low signal to noise ratio. 
The lack of a ^-dependent acceptance correction can also have an impact on the 
normalization in 3-particle correlations. To improve this correction we use a single 
particle dependent acceptance correction. This correction is done separately for each 
charge and magnetic field setting. Figure 13.2.21 shows examples of the acceptance 
for two polarities of magnetic field settings for our associated particles in central 
Au+Au collisions. The acceptance correction plots are binned in very narrow 
bins to make any binning effect negligible. Particles are weighted by the reciprocal 
of the dependent acceptance. The dependent acceptance is normalized such that 
the average is one, since we have already corrected for the number of particles with 
the overall 0-averaged detector efficiency. Because high particles curve little in 
the STAR magnetic field, the trigger particle acceptance goes to zero at the TPC 
sector boundaries. To avoid large corrections we do not use trigger particles that 
have less than a 10% acceptance. Since the quantities are per trigger this has no 
effect on our overall efficiencies. The number of trigger particles is accumulated by 
the 0-dependent correction factor, and is used in the final normalization of the per 
trigger correlation functions. 

3.3 Background Construction 

We are using charged particles in the STAR TPC. Correlations are between a 
trigger particle of 3 < < 4 GeV/c and two associated particle of 1 < < 2 
GeV/c. The data used in this thesis are from 2 million pp events from the second 
year of RHIC running, 6.5 million d+Au events used from third year of RHIC running. 
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Phi effiency correction 1.00Pt2.00 FF Positive 5-0% inciusive Phi effiency correction 1.00Pt2.00 RFF Positive 5-0% Inciusive 




Figure 3.3. TPC acceptance in for particles of 1 < < 2 GeV/c 
with ibO.5 Tesla magnetic field. Positive particles are shown in black and 
negative particles are shown in red. Plots are for 0-5% Au+Au collisions. 
Left: Magnetic field with postive polarity. Right: Magnetic field with 
negative polarity. Errors are statistical. Plots for other centralities and 
for 3 < pt < 4: GeV/c particles can be found in the appendix. 



and 12 million minimum bias and 19 million central triggered events from the fourth 
year of RHIC running. The number of events listed are good events after vertex 
cuts. The tracks used all have a distance of closest approach (DCA) to the primary 
vertex less than 2 cm. They have at least 20 points fit by the track reconstruction 
software and at least 51% of the maximum possible number of fit points (determined 
by the detector geometry for the particular track) to avoid split tracks. All tracks are 
within \ri\ < 1.0. The events were divided in centrality bins based on the so-called 
reference multiplicity. In Au+Au collisions, the reference multiplicity is the number 
of tracks with DCA < 1 cm, number of fit points greater than 10, and \ri\ < 0.5 
(the STAR reference multiplicity). In d+Au collisions, the reference multiplicity 
is the number of tracks with DCA < 2, number of fit points great than 15 and 
|?7| < 1.0. The centrality divisions are shown in Table [3?T1 The analysis was performed 
in these individual centrality bins and then combined, weighted by the number of 
trigger particles, into larger centrality bins. This was done to reduce non-Possion and 
multiplicity dependent effects that increase with the size of the centrality bin. 
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Table 3.1 

Centrality bins definitions used in this thesis in d+Au and Au+Au col- 
hsions. The minimum and maximum are apphed on the reference multi- 
phcity. See text for the definition of reference multiphcity in d+Au and 
Au+Au colhsions. 



Centrahty 


Minimum 


Maximum 


d+Au 100-20% 





17 


d+Au 20-10% 


18 


26 


d+Au 0-10% 


27 


100 


Au+Au 70-80% 


15 


31 


Au+Au 60-70% 


32 


57 


Au+Au 50-60% 


58 


96 


Au+Au 40-50% 


97 


150 


Au+Au 30-40% 


151 


222 


Au+Au 20-30% 


223 


319 


Au+Au 10-20% 


320 


441 


Au+Au 5-10% 


442 


520 


Au+Au 0-5% 


521 


1000 



Figure [33] shows an example raw 3-particle correlation function, J3(A0i, A</)2), 
where A(/)j = (pi — (pxrig- The raw 3-particle correlation function contains our desired 
jet-like 3-particle correlation and background terms. These background terms are 
discussed in detail below. 

3.3.1 Hard-Soft Background 



The two-particle correlation function is given by (as discussed in Chapter 3), 

J2(A0) = J2(A0) - a5,„eF2(A0) (3.14) 
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Figure 3.4. Raw 3-particle correlation function for 0-12% central Au+Au 
collisions. 



where 

F2'\A(j)) = 1 + 2v^v^^cos{2A(j)) + 2<t;^cos(4A0) (3.15) 

is the 2-particle flow modulation up to the fourth harmonic and Bine = dNinc/d{A(j)). 
Here our background event is an inclusive event, a minimum bias event of the same 
centrality. The normalization factor a is used to adjust the level of the background due 
to trigger bias effects. In the 3-particle correlation analysis a will not be normalized 
to 2-particle ZYAl as is done in the 2-particle analysis (presented in Chapt. 3). Here 
we will use the assumption the 3-particle jet-like correlation is positive definite and 
use 3-particle ZYAM (zero yield at minimum). This can provide a better constraint 
on the 2-particle background since we have additional information in the 3-particle 
correlations. The systematics on this assumption will later be discussed. An example 
raw 2-particle correlation J2{A(j)) is shown in figures [375] (left) . in filled symbols. The 
normalized background from mixed events is shown by the solid black line with flow, 
aBincF2{A(f)). The mixed events were constructed by mixing the trigger particle with 
inclusive events. The flow, both f 2 and v^, is added pairwise during the event mixing. 
The V2 is the average of the reaction plane and 4-particle cumulant measurements. 
The f4 is taken to be 1.15f|. The minipanel shows the background subtracted 2- 
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particle correlations, J2. The right panel shows the hard-soft background term. This 
term is given in equation 13.41 This is a folding of the background subtracted 2- 
particle correlation and its background. This background term removes 2-particle 
correlations, where one particle is correlated with the trigger and the other is only 
correlated via the reaction plane, from our 3-particle correlation. In real data analysis 
this term becomes, 

J2 ® ^2 = J2(A0i)a5,„e[l + Ft^^''\^(t>2)] + J2(A02)a5,„e[l + Ft^^^\^<Pi)] 
= [J2(A0i) - aB,^,Fr-'''^'\A<p^)]aB,nc[l + Ft''^'\A<P2)] 
+ [J2(A02) - aB,^,F^"''^'\AcP,)]aB,^,[l + F^'''^\A<P,)]. 

(3.16) 

Then we fold the jet-like 2-particle correlation with the flow modulated background. 
To obtain f 2 we make the assumption that the trigger-associated jet-like correlated 
pair are correlated with the reaction plane with the same flow as the trigger particle 
(^2^* = v^^^^). The systematics on this assumption will be discussed later. 



3.3.2 Soft-Soft Background 

Figure [3^ shows the soft-soft background term. This term contains all correlations 
between the two softer particles that are independent of the trigger particle. In pp and 
d-|-Au collisions, this term is composed of minijets, background jets, and decays. In 
Au-|-Au collisions, there is also a contribution from the flow correlation between the 
two softer particles. This flow component is the dominant part of the soft-soft term 
in Au-|-Au collisions. The soft-soft term is constructed by mixing a trigger particle 
from one event with pairs of softer particles from another (inclusive) event. It is the 
inclusive event particle pair density with respect to a random trigger, 

^--lA^^Afe)^ "^^;^;;^-;)' (3.17, 

The flow component already contained in this term is, 

jr'^'°-(A0i,A02) = - ^) [1 + 2v^i\^? cos2(A0, - A<j>,) 
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Figure 3.5. Left: Raw 2-particle correlation function is shown in red, 
J2. Normalized background from mixed events with v2 and f4 added, 
aBincF2, is shown in black. Normalization is done by 3-particle ZYAM. 
Background subtracted 2-particle correlation is shown in the minipanel, 
J2. Right: Hard-soft background term, J2 ® -82- Plots are from 0-12% 
central Au-|-Au collisions at ^/sNlv = 200 GeV/c. Errors are statistical. 



+24^^f4^^ cos4(A0i - A02)] (3.18) 

and is shown in Fig. 13.61 (right). This was constructed by mixing the trigger particle 
with two different inclusive events and adding the flow in tripletwise. There is a 
{Ninc{Ninc — 1))/ {N-HN-^I) factor applied so that the number of associated pairs is 
the same as if they came from the same event. The soft-soft background is scaled 
by a'^b where a is the scaling factor used in the 2-particle correlations and accounts 
for triggering effects, as discussed in Chapt. 3, and b accounts for non-Poisson ef- 
fects. The factor a is determined by a 3-particle ZYAM assumption. The factor b is 
required because the event multiplicity distributions are not Poisson. The factor b is 
determined by the ratio of the deviation from Poisson of the triggered events to the 
inclusive event. 



61 



This method of obtaining b makes the assumption that the underlying background 
and the triggered events have the same level of deviation from Poisson. 




A(b =(b -(b A(b =(b -(b 

^1 ^1 ^Trig ^1 ^1 ^Trig 

Figure 3.6. Left: Soft-soft background term, o?hJl^^. Right: Flow compo- 
nent of the soft-soft background term, a?hJ^^'^'^^°^ . Plots are from 0-12% 
central Au+Au collisions at ^/snn = 200 GeV/c. 



3.3.3 Trigger Flow Correlation 

There are additional flow correlations that are not preserved in the event mixing 
and thus must be subtracted by hand. There is the flow correlation between the trigger 
particle and the two associated particles, where both of the associated particles come 
from the underlying background in the triggered event. The total trigger flow is, 

Jr = ^^^^%^i^3(A0i, (3.20) 

with 

F3(A0i,A02) = F2(A0i) + F2(A02) 

+2i;*"^4'^t;f cos (2 A01 - 4A02) 
+2vi"'v^\^2^cos{AA(j)i - 2A02) 

+2i;f ^4'^ cos (2 A01 + 2A02). (3.21) 
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Figure 13.71 (left) shows the flow correlation with the trigger particle up to order f 2. 
The central panel shows the flow correlation with the trigger particle of order f| 
where assumed proportional to The sum of these two panels is Jj^^{A(f)i, A02)- 
These terms are constructed by mixing the trigger particle with associated particles 
from two different inclusive events. The flow terms are added in tripletwise with 
the V2 from the average of the reaction plane and 4-particle cumulant measurements 
and V4 = 1.15^2. Since the associated particles come from two different inclusive 
events there is a factor of {Ninc{Ninc — l))/(^inc^jnc) applied so that the average 
number of associated pairs is same as the average number of associated pairs when 
both associated particle come from the same event. As with the soft-soft term, the 
trigger flow terms are scaled by a'^b. 



Ad) =(1) -0 

1 1 ^Trig 



• f 





Ad) =(1) -6 
1 1 ^Trig 



Figure 3.7. Left: Trigger flow correlation up to order f|. Center: Trigger 
flow correlation of order Right: Trigger flow correlation with non- 
flow on soft-soft, Jj^ ■ Plots are from 0-12% central Au-|-Au collisions at 
= 200 GeV/c. 



This would be exact if the particles were only correlated by flow. However, there 
is an additional term required due to the flow correlation of the non-flow structure 
on the soft-soft term (such as minijets, decays and background jets) with the trigger 
particle. We need to correct for this effect in the trigger flow. To first order this can 
be approximated as, 

^CF \ Jr(A0i,A02 



jinc,flow 

Jr, 



(A0i,A02) 



F3(A0i,A02) 



(3.22) 
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An example of this term is shown in Fig. 13.71 (right). Since the ratio of the soft-soft 
term and the flow between the two soft particles is close to 1, this term is small and 
the first order correction is sufficient. 



3.3.4 Background Subtraction 



An example background subtracted 3-particle correlation is shown in Fig. l3.8[ left. 
This is obtained from the raw signal through the subtraction of the background terms 
discussed by, 

MA(j)i, Acj)2) = J3(A0i,A02) 

-aJ2(A0i)5,„eF/^*'(')(A02) - aUA<P2)B,^,F^''''~'\A<f),) 
-a%[Jt%A<f)i, A<P2) 

-{F3(A0i, A02) + F^'^^XAcP,, A02)}] 



(2vr)^ 

(3.23) 

where a is determined such that J3(A0i, A02) is ZYAM and . To obtain the ZYAM, 
the lowest 10% of the bins are used and a is adjusted till these bins average zero. The 
bins used for the minimum are recalculaed for each values of a. The factor b is there 
because the event multiplicity distributions are not Poisson. 

One item of note in the background subtraction is the partial cancellation of the 
flow. The flow term of f J"^f2*^ and v'^"^v'^\ where i = 1,2, partially cancel between 
the flow subtracted to form 2-particle jet-like correlation used to construct the hard- 
soft term and the terms in the trigger flow. If the events are Poisson, then these 
terms will entirely cancel. This partial cancellation gives a much smaller uncertainty 
on the sum of the hard-soft term and the trigger flow terms (Fig. 13.81 right) due to the 
uncertainty in the elliptic flow measurement than they have individually. This allows 
us obtain a significant background subtracted 3-particle signal even in bins where the 
uncertainties on the hard-soft and trigger flow terms due to the uncertainty on the 
f 2 measurement are larger than our signal. 
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are from 0-12% central Au+Au collisions at ^/sNN = 200 GeV/c. 



3.4 Results 

3.4.1 Correlation Functions 

Figure 13.91 shows the background subtracted 3-particle correlations for pp, d+Au, 
four centralities of Au+Au from minimum bias data and ZDC triggered Au+Au. The 
pp and d+Au background subtracted 3-particle correlations are very similar. Both 
show four peaks. A peak is at (0, 0) from instances where both associated particle are 
on the near-side. Another peak is at (tt, tt) from instances where both particles are 
on the away-side. This peak is slightly elongated along the diagonal. This elongation 
is qualitatively consistent with broadening which is due to nonzero initial total 
transverse momentum of the colliding partons. This results in di-jets that are not 
exactly back-to-back. This effect is enhanced by selecting high pt trigger particles 
because the selection preferentially picks up those parton scatterings that have a large 
total transverse momentum. The fc^ is along the diagonal because it displaces the 
away-side from tt but all the jet particles on the away-side stay close together. The 
peaks at (0, tt) and (tt, 0) are also elongated for the same reason. 
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In the peripheral Au+Au colhsions, additional broadening is seen along the diag- 
onal. This additional elongation could be due to jets deflected by radial flow, path 
length dependent energy loss, or large angle gluon radiation. In the 30-50% Au+Au 
centrality bin, potential peaks start to develop on the off-diagonal, suggesting a small 
contribution from conical emission. The on-diagonal elongation continues into this 
bin. There is also on-diagonal elongation of the near-side. This is because we are 
using the trigger particle azimuthal angle as a proxy for the azimuthal angle of the 
jet-axis. When the trigger paritcle is to one side of the jet-axis, we are more likly to 
get both associated particle to the other side of the jet-axis. This also has the effect of 
slanting near-away peaks. Figure [3. 101 is a cartoon to help visualize this effect. In the 
10-30% Au+Au bin, definite off-diagonal peaks are visible. The 0-10% Au+Au bin 
seems to have a simialar signal; however the statistics are very poor. For increased 
statistics in central collisions we have taken enriched central data samples afforded by 
an online ZDC trigger. This trigger provides a 0-12% most central Au+Au collsions 
with an order of magnitude more statistics. With this data sample, very distinct 
off-diagonal peaks can be seen, providing unambiguous evidence for conical emission. 

To examine the results quantitatively, we calculated the average yield, of the 
background subtracted 3-particle correlation, in boxes centered at regions of interest. 
Figure 13.111 (left) shows a cartoon with the locations of the boxes. The right panel 
shows the average yields at these locations, the near-side (0,0), the away-side {n,7i), 
on-diagonal (vr ± 1.42, tt ± 1.42) and off-diagonal (vr ± 1.42, tt + 1.42) plotted as a 
function of centrality. The near-side yield increases with centrality. The away-side 
yield increases through pp, d+Au and peripheral Au+Au collisions but seems to level 
out in mid-central to central Au+Au collisions. The on-diagonal and off-diagonal 
yields are consistent with zero in pp, d+Au, and peripheral Au+Au collisions. In 
mid-central Au+Au collisions, they are significantly above zero. The on-diagonal 
peak is consistently larger than the off-diagonal peak. 

We obtain more quantitative information by studying the projections of the back- 
ground subtracted 3-particle correlations. We expect two side peaks in the off- 
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Figure 3.9. Background subtracted 3-particle correlations in (from right to 
left, top to bottom) pp, d+Au, Au+Au 50-80%, Au+Au 30-50%, Au+Au 
10-30%, Au+Au 0-10% and ZDC triggered Au+Au 0-12%. collisions at 
= 200 GeV/c 



diagonal projection for conical emission. Figure 13.131 shows on-diagonal and off- 
diagonal projections of the background subtracted 3-particle correlations shown in 
Fig 13. 9[ These projections are done by projecting strips of width 0.7 radians about 
the diagonal (off-diagonal) to the diagonal (off-diagonal) for particles on the away- 
side, a diagram showing the projected regions is shown in Fig. /reffig:ProjDiagram. 
Since we have finite binning and we are not projecting along one of the axes, bin width 
effects will appear in these projections if we do not take care of them. This is removed 
by sampling each individual (A0i, A02) bin randomly 100 times (i.e. 1/100*^ of the 



Figure 3.10. Useful for visualizing the on-diagonal elongation and the 
tilting of the near-away peaks. The black line is used to represent the 
jet-axis. The red line is used to represent the trigger particle. The bumps 
represent the near-side and away-side distribution. 




Figure 3.11. Left: Cartoon with approximate positions of squares used 
for the yields. Right: Average yield in 0.7x0.7 squares centered on the 
near-side (0, 0), the away-side (tt, tt), on-diagonal (vr ± 1.42, tt ± 1.42), and 
off-diagonal (tt ± 1.42, tt =p 1.42). Npart is the number of participants. 
The solid errors are statistical and the shaded are systematic. The ZDC 
triggered 0-12% most central Au+Au coUions point is shifted to the left 
for clarity. Collisions are at ^/sNN = 200 GeV/c. 
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bin content is projected at 100 random positions within each bin). This removes the 
jagged effects due to finite bin width but does leave us with bin smearing of up to 
half a bin width. 

5 

-&* 
O 

4 
3 
2 



A*, 

Figure 3.12. Regions that are projected for the on-diagonal and off- 
diagonal projections are shown in blue and red respectively. 

As seen from Figure 13. 131 in pp and d+Au collisions the on-diagonal (blue) and off- 
diagonal (red with yellow systematic error boxes) projections look very similar, both 
are single peaked. The only difference is the on-diagonal projection is wider. This 
is likely due to kx broadening as previously discussed. The most peripheral Au-|-Au 
bin (50-80%) has a very broad on-diagonal peak. The off-diagonal projection shows a 
central peak and two symmetric side peakj^. These side-peaks are not very significant, 
perhaps indicating only a small contribution from conical emission to this centrality. 
The next centrality bin (30-50% most) shows an on-diagonal peak that is very broad 
and consistent with being flat on top. The off-diagonal projection shows side peaks 
that are larger than the 50-80% centrality bin, however the systematic errors are 
also larger and the peaks are only about two sigma effect. In the 10-30% Au+Au 
centrality bin, the on-diagonal projection has become dipped at [A0i -|- A02]/2 = vr. 
The off-diagonal projection shows significant side peaks, evidence for conical emission. 

"'^The 3-particle correlation functions are symmetric with respect to A<pi and A<p2 by construction so 
the off-diagonal projections are always symmetric, except for effects from the randomization and the 
bins are not set such that tt is a division. The on-diagonal projection has no such forced symmetry. 
The statistical errors in the projections do not have a correction for the double counting of associated 
pairs. 
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The 0-10% most central Au+Au collisons look very similar although the statistical 
errors are much larger. The 0-12% ZDC triggered central Au+Au data again looks 
similar (but with smaller errors) with a dip in the on-diagonal projection, a centeral 
peak and two symmetric side peaks in the off-diagonal projection. These side peaks 
are consistent with conical emisson in the central Au+Au data. 
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Figure 3.13. Projections of strips of full width 0.7 radians on the away- 
side. On-diagonal projections (blue) and off-diagonal projection (red) 
in centralities (from left to right, top to bottom) 'p'p^ d+Au, Au+Au 50- 
80%, Au+Au 30-50%, Au+Au 10-30%, Au+Au 0-10% and ZDC triggered 
Au+Au 0-12% collisions at ^/sNN = 200 GeV/c. Error bars are statistical 
errors. Systematic errors on the off-diagonal projections are shown in the 
yellow band. Dashed black lines are at zero. The on-diagonal projections 
with systematic errors can be found in the appendix. 



There is extra signal on the on-diagonal projections when compared to the off- 
diagonal projections. This is shown by the difference in the projections in Fig. 13.141 
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The difference plotted is the on-diagonal projection minus the side Gaussians obtained 
from the fit to the off-diagonal projection. In pp and d+Au, there is little signal in 
the off-diagonal projections in the region of the side peaks so this difference is very 
similar to the on-diagonal projection. In central and mid-central Au+Au collisions, 
there is a central wide peak that is rather ffat. There is also some small contribution 
on the outside edges. The peaks seen in the difference could be due to many differ- 
ent physics mechanism such as: large angle gluon radiation, jets deffected by radial 
ffow or preferential selection of particles due to path-length dependent energy loss, 
hydrodynamic conical ffow generated by Mach-cone shock waves that couple with 
flow [69], and/or jets for which the away-side jet undergoes relatively little medium 
modification. 

3.4.2 Emission Angle and Associated Pt Dependence 

We can extract the conical emission angle by fitting the off-diagonal projections to 
a central Gaussian and symmetric side Gaussians. Table [321 gives the angles obtained 
from the fits for mid-central and central Au+Au collisions. They are also plotted in 
figure 13.151 If the emission angle is the Mach-cone angle, we may use this angle to 
extract the average speed of sound of the medium, Cg, from the obtained Mach-cone 
angle 9^^ , by. 



This is the average speed of sound, averaging over the entire time evolution of heavy 
ion collisions, which may undergo several stages: QGP, phase transition, and hadronic 
stages. Expected sound velocities for different phases are listed in table 13.31 This 
equation gives Cg = 0.15c for 6'^^ = 1.42 radians. One should be cautious about this 
number because it neglects effects of hydrodynamics (other than Mach-cone) and 
expansion. The relation between the emission angle in data and the speed of sound 
in the medium could be very complicated and is a subject of on-going theoretical 
work. 




(3.24) 
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m^-n)±{A^^-n))l2 ((A(|)^-7t)±(A(|)^-7i))/2 ((A(^^-7i)±(A(l)^-7i))/2 



Figure 3.14. Difference between in the on-diagonal projection and the 
side Gaussians from fit to the off-diagonal projection. Projections are 
in strips of full width 0.7 radians on the away-side in centralities (from 
left to right, top to bottom) pp, d+Au, Au+Au 50-80%, Au+Au 30- 
50%, Au+Au 10-30%, Au+Au 0-10% and ZDC triggered Au+Au 0-12% 
collisions at ^/sNN = 200 GeV/c. Error bars as statistical from the on- 
diagonal projection only. Systematic errors are shown in the yellow band. 
Dashed black lines are at zero. 



It is important to investigate the emission angle as a function of associated particle 
Pt- For conical emission from a Mach-cone the angle is expect to remain constant as a 
function of associated particle px, since the speed of sound is an intrinsic property of 
the medium and independent of the partonic momentum. Simple Cerenkov radiation 
models predict a sharply decreasing angle as a function of associated particle pr [71] . 
This Pt dependence can be used to distinguish conical emission from Mach-cone and 
Cerenkov radiation. 
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Table 3.2 

Angles from Gaussian fits to the off-diagonal projections. 



Centrality 


Angle 


Statistical Error 


Systematic Error 


AuAu 30-50% 


1.43 


±0.01 


+0.04 
-0.03 


AuAu 10-30% 


1.42 


±0.01 


+0.04 
-0.03 


AuAu 0-10% 


1.40 


±0.03 


+0.07 
-0.05 


AuAu ZDC 0-12% 


1.39 


±0.01 


+0.04 
-0.04 



I 1.9 

CO 



< Au+Au 0-12% (shifted) 
I Au+Au 30-50%, 10-30% and 0-10% 

l^-s^ 1.4210.02 
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Figure 3.15. Angles from Gaussian fits to the off-diagonal projections. The 
statistical errors (solid) are the fit errors. Systematic errors are shaded. 
The blue line is a fit of the points to a constant, yielding 1.42 ± 0.02 
(fit error using quadrature sum of statistical and systematic errors on the 
points). Npart is the number of participants. The ZDC triggered 0-12% 
most central Au±Au point is shifted to the left for clarity. The dashed 
line is at n/2. The collisions are at y^sJjN = 200 GeV/c. 



Figure [3. 161 shows background subtracted 3-particle correlations in ZDC triggered 
0-12% Au±Au collisions for different associated particle pt bins. The pt dependent 
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Table 3.3 

Expected speeds of sound in different phases of matter [781179] . 



Pfiase 


Speed of Sound 


QGP 




Mixed Pliase 





Resonance Gas 


0.47c 




Figure 3.16. Background subtracted 3-particle correlations for ZDC trig- 
gered 0-12% most central Au+Au collisions at ^/sNN = 200 GeV/c. Trig- 
ger particle is 3 < pt < 4 GeV/c. Associated particle transverse momen- 
tum from left to right top to bottom is 0.5 < pr < 0.75, 0.75 < pr < 1, 
1 < Pt < 1.5, 1.5 <pt <2, and 2 < px < 3 GeV/c. 

3-particle correlation plots were not corrected for acceptance on a seperate basis for 
the different magnetic field settings. However, the results for the 1 < px < 2 GeV/c 
associated particles have very little dependence on whether or not this correction 
was done for both field setting together or seperately for the different fields. The 
magnitude of the off-diagonal peaks decreases with increasing associated particle pt as 



74 



expected. Figure [3TT71 shows the on-diagonal and off-diagonal projections of the plots 
in Fig. I3.16[ The on-diagonal projections with the conical emission from the fit to the 
off-diagonal projection removed are shown in Fig. 13.181 From Figures 13.171 and 13.181 
it can be seen that the relative difference between the on-diagonal signal and the off- 
diagonal signal increases with associated particle pt- Figure [3. 191 shows the angle from 
the Gaussian fits to the off-diagonal projections as a function of associated particle pt- 
The angle is consistent with remaining constant as a function of associated particle pt 
as expectd for Mach-cone emission. It is inconsistent with a sharply decreasing angle 
as a function of associated particle px as predicated by simple Cerenkov radiation [71j. 
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Figure 3.17. Away-side projections of background subtracted 3-particle 
correlations for ZDC triggered 0-12% most central Au-|-Au collisions at 
y/SNN = 200 GeV/c. The on-diagonal projection is in blue and off- 
diagonal projection in red. Yellow bands represent the systematic error 
on the off-diagonal projection. Trigger particle is 3 < p-r < 4 GeV/c. 
Associated particle transverse momentum from left to right top to bot- 
tom is 0.5 < Pt < 0.75, 0.75 < pt < I, I < Pt < 1-5, 1.5 < pt < 2, 
and 2 < Pt < 3 GeV/c The projections are of strips of full width of 0.7 
radians. 
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m-m^^%-^m m-m^<^^-^)V2 m^-Tzm^-^w 



Figure 3.18. On-diagonal projections with the side Gaussians from the 
off-diagonal projection removed. The red curves are the Gaussians that 
were subtracted. Yellow bands represent the systematic error on the off- 
diagonal projection. Trigger particle is 3 < < 4 GeV/c. Associated 
particle transverse momentum from left to right top to bottom is 0.5 < 
Pt < 0.75, 0.75 < pt < I, I < Pt < 1-5, 1.5 < px < 2, and 2 < < 3 
GeV/c The projections are of strips of full width of 0.7 radians. Plots are 
for 0-12% ZDC triggered Au+Au collisions at y/sj^ = 200 GeV/c. 



3.4.3 Widths 

The widths of the conical emission peaks are shown as a function of centrality 
in Fig. I3.20[ top left. The widths are constant with being independent of centrality 
or increasing with centrality. An increase in the width with centrality could point 
towards more dispersion of the cone through more interactions with the medium in 
central collisions. The top right and bottom left plots show the widths of the central 
away-side peak in the off-diagonal and on-diagonal projections, respectively. There 
is not significant centrality dependence of central peak in the off-diagonal projection. 
The width of the central peak in the on-diagonal projections is the width of the peak 
shown in Fig. I3.14[ This peak appears to increase in width as the centrality increases. 
The interpretation of the broadening is complicated as previously discussed and could 
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0-12% 1.36±0.03 



I I I I I 

2 5 3 
P,(GeV/c) 

Figure 3.19. Angles from fits to off-diagonal projections as a function 
of associated particle pt for ZDC triggered 0-12% most central Au+Au 
collisions at ^/sJfJv = 200 GeV/c. Trigger particle is 3 < < 4 GeV/c. 
The solid errors are statistical and the shaded are systematic. The blue 
line is a fit to a constant, yielding yielding 1.36 ib 0.03 (fit error using 
quadrature sum of statistical and systematic errors on the points). 

come from many different physics mechanisms. The bottom right plot shows the 
difference in the widths of the on-diagonal projection and the off-diagonal projection. 
This difference increases with centrality. Again the interpretation is complicated but 
the broadening of the on-diagonal central peak relative to the off-diagonal central 
peak could be due to interactions between away-side jet and the medium. 

Fig. 13.211 shows the centrality dependence of widths of the near-side peak, in both 
the off-diagonal and on-diagonal projection, top left and right, respectively. The 
near-side peak width is consistent with being independent of centrality. This is not 
unexpected if most of the trigger particles come from surface emission. The on- 
diagonal width is wider than the off-diagonal width because we are using the trigger 
particle as a proxy for the jet-axis. If the trigger particle is to one side of the jet- 
axis than both associated particles are more likely to be to one side of the trigger 
particle (which puts them on-diagonal) and more likely to be further from the trigger 
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Figure 3.20. Centrality dependence of the width (a of Gaussians) from fits 
to away-side projections of the background subtracted 3-particle correla- 
tions. Top Left: Width of the side Gaussians in the off-diagonal projec- 
tion. Top Right: Width of the central peak in the off-diagonal projection. 
Bottom Left: Width of the central peak in the on-diagonal projection 
after subtraction of the side Gaussians from the fit to the off-diagonal 
projection. Bottom Right: Difference in the widths of the central peaks 
done as sJ(tI„_ diagonal - <^off -diagonal- The uumbers indicate the constant 
fit results. The solid errors are statistical and the shaded are systematic. 



particle than the jet-axis (which results in the broadening). The difference in the 
widths is shown in the bottom panel. This difference is consistent with no centrality 
dependence. 

Figure 13.221 shows the widths of the away-side peaks as a function of associated 
particle px- Although the errors are large, the conical emission peaks (left) seem to 
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Figure 3.21. Centrality dependence of the width (cr of Gaussians) from 
fits to near-side projections of the background subtracted 3-particle 
correlations. Top Left: Off-diagonal peak width. Top Right: On- 
diagonal peak width. Bottom: Difference in the peak widths as 
^(^In-diagonal " (^If f -diagonal- The uumbers indicate the constant fit re- 
sults. The solid errors are statistical and the shaded are systematic. 



have associated particle px independent widths. The width of the central peak (right) 
is consistent with decreasing with associated px- 

Figure 13.231 shows the widths of the near-side peak as a function of associated par- 
ticle pt- The widths of the on-diagonal projection appear to decrease with associated 
particle pt- This is probably because the high pt associated particles will be more 
aligned with the jet-axij^. The off-diagonal near-side projection widths are consis- 

^This is because the average momentum perpendicular to the jet-axis (jt) is independent of the 
associated particle pt- For higher particles this is a smaller fraction of their transverse momentum 
so it leads to a smaller angular deviation from the jet-axis. 
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Figure 3.22. Width (a of Gaussians) of away-side peaks from the off- 
diagonal projections of the background subtracted 3-particle correlations 
as a function of associated particle pr- Left: Width of the side Gaussians. 
Right: Width of the central Gaussian. Plots are for 0-12% ZDC triggered 
Au+Au at y^s]^ = 200 GeV/c with 3 < p^"^ < 4 GeV/c. The numbers 
indicate the constant fit results. The solid errors are statistical and the 
shaded are systematic. 



tent with either a decrease with associated particle p-p or independent of associated 
particle pt- There is no significant change in the difference of the widths. 

3.4.4 Trigger Pt Dependence 

We can also look at how the 3-particle correlation changes with trigger particle pt- 
Figure [3.241 (left) shows the background subtracted 3-particle correlation for triggers 
of 4 < Pt < Q GeV/c and Q < pt < 10 GeV/c. The result looks very similar to that 
for 3 < < 4 GeV/c. There does appear to be a little less on-diagonal broadening 
on both the near-side and the away-side for the higher pt trigger. This is expected 
because the azimuthal angle of the higher pt trigger is a better proxy for the jet-axis. 
Figure 13.241 also shows the on-diagonal and off-diagonal projections. The yields and 
the angle of the side peaks in the off-diagonal projection are shown in Fig. 13. 251 The 
yields increase with increasing trigger particle pt- The away-side yield increases more 
quickly then the off-diagonal peak yield. This is consistent with more of the away-side 
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Figure 3.23. Width (cr of Gaussians) of near-side peaks from projections 
of the background subtracted 3-particle correlations as a function of as- 
sociated particle pr- Top Left: Off-diagonal peak width. Top Right: 
On-diagonal peak width. Bottom: Difference in the peak widths as 
\l (^In-diagonal " (^If f -diagonal- Widths are from Au+Au 0-12% ZDC trig- 
gered data at ^/snn = 200 GeV/c. The numbers indicate the constant fit 
results. The solid errors are statistical and the shaded are systematic. 



jet punching though the center for higher energy jets. The angle of the side peaks in 
the off-diagonal projection is consistent with no trigger particle pr dependence. 

3.5 Systematic Uncertainties 

Systematic error bars have been shown in different plots in the results section. 
This section discusses what goes into the systematic errors. The systematics have 
been rigorously studied. There are two dominant sources of systematic uncertainties. 
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Figure 3.24. Three-particle correlations with trigger particles of 4 < 
p^*^ < 6 (top) and 6 < p^^^^ < 10 (bottom) with associated particles 
of 1 < < 2 GeV/c for 0-12% ZDC triggered Au-hAu collisions at 
y/SNN = 200 GeV/c. Left: Background subtracted 3-particle correla- 
tion. Right: On-diagonal projection (blue) and off-diagonal (red). Error 
bars are the statistical error. Yellow band is the systematic error on the 
off-diagonal projection. 



These are normalization and uncertainty on the flow measurement. Other sources of 
systematic error include V2 of the jet, effect of flow fluctuations, the parameterization 
of f4, and finite multiplicity bin width effects. There is also an overall 10% systematic 
uncertianity from the uncertianity in the efficiency. In the following sections, we 
discuss the sources of systematic uncertainty in detail. 
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Figure 3.25. Left: Average yield in 0.7x0.7 squares centered on the near- 
side (0,0), the away-side (vr,7r), on-diagonal (vr ± 1.42, tt ± 1.42), and 
off-diagonal (vr ± 1.42, vr =l= 1.42) as a function of trigger particle Pt- Right: 
Angles from fits to off-diagonal projections. The solid errors are statistical 
and the shaded are systematic. Plots are for 0-12% ZDC triggered Au-|-Au 
collisions at ^/snn = 200 GeV/c. 



3.5.1 Background Normalization 

We have two factors that are used in our background normalization. Both of 
these parameters have been explored in parameter space. The default values have 
been chosen based on our best knowledge of jet-like correlations at RHIC 

The normalization factor a was at first chosen to be the same as is used in 2- 
particle correlations. The assumption used was that the 2-particle correlated yield 
was zero at A0 = 1 (the minimum of the background subtracted 2-particle correlation 
in Au-|-Au collisions falls at 1). This was always known as a lower limit on the yield 
and an upper limit on a but in 2-particle correlations there is not enough information 
to do any better. This is an upper limit because our signal represents the number of 
particles (in a particular range) associated with a high pjn trigger particle which is 
positive definite. This can now be improved upon in 3-particle correlations. The 3- 
particle correlation signal should be positive definite for the same reason, the number 
of pairs associated with a high px trigger particle should be positive definite. We 
therefore determine a using the assumption that the 3-particle correlation signal is 
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zero yield at minimum. This is done by selecting the 10% lowest bins (58 of 24 x 24) 
and adjusting a until the average content of these bins is zero. The lowest bins were 
redetermined for each adjusted value of a. For the systematic error assigned to our 
normalization factor a, we use our known upper limit, the value from 2-particle ZYAl. 
The difference between the default value and the value obtained from the 2-particle 
ZYAl is considered to be our uncertainty. Therefore for our lower limit we subtract 
this difference from the default value. Figure 13.261 shows the background subtracted 

2- particle correlations when a is obtained from 2-particle ZYAl (our upper systematic 
limit) and from when a is twice the default value minus the 2-particle ZYAl value 
(our lower systematic limit) in the left and right panels, respectively. Table [3^ shows 
the default values of a for each centrality bin along with the systematic errors applied. 

In the beginning of this analysis, the normalization factor b was not applied, the 
equivalent of 6 = 1. It was obvious in d-l-Au collisions that this factor was neces- 
sary, for without this factor there was a large pedestal in the background subtracted 

3- particle correlation signal. This factor is necessary because the multiplicity distri- 
butions are non-Poisson. If the events were Poisson than {N{N — 1)) = (A^)^. From 
data we know this is not true. We did not require the events to be Poisson in our 
analysis, but assumed that {N{N — 1))/{N)'^ is the same for the triggered events, 
inclusive events and the underlying background. This is a less stringent requirement. 
With this assumption the soft-soft and trigger flow terms are scaled by a^. However, 
we have found that {N{N — 1))/ (A^)^ is not the same for the triggered events and the 
inclusive events so we introduce 6 as a correction, where 

{NtrigjNtrig-'i-)) 

b= • (3.25) 

We are therefore assuming that the underlying background deviates from Poisson 
similar to the deviation in triggered events. 

The uncertainty applied to b comes from the uncertainty in the 3-particle ZYAM. 
This takes care of our uncertainty on ZYAM due to the number of bins used and 
gives us a relative change between a and b. We use the default value of a that was 




Figure 3.26. Background subtracted 3-particle correlations for the un- 
certainty on normalization factor a. Left: Upper limit on a where it is 
obtained from 2-particle ZYAl. Center: Lower limit on a where it is the 
default from 3-particle ZYAM minus the difference between the upper 
limit and the default. Right: Off-diagonal projection from default with 
systematic uncertainty from the uncertainty on a shown in shaded band. 
From top to bottom plots are pp, d-l-Au, Au-|-Au 50-80%, Au-|-Au 30- 
50%, Au+Au 10-30%, Au+Au 0-10%, and ZDC triggered Au+Au 0-12% 
collisions at ^Jsnn = 200 GeV/c. 
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Table 3.4 

Values for normalization factors a and h with systematic errors. ZDC 
triggered Au+Au collisions are divided into centrality bins as if they were 
minimum bias events. The cross section numbers in quotes correspond to 
the cross sections as if events were minimum bias but not the actual cross 
sections for these events. 



Collision 


Centrality 




h 




0-100% 


687 ± 1 040 

\J • \JK-3 f — 1 — -I- • \J^I.\J 


0.81 + 0.60-0.82 




20-100% 


1 001 ± 039 


1 9fi4 4- n m 


n (199 




10-20% 


652 ± 177 


1 QzL9 -1- n (T^n 


— n nzLQ 


d+Au 


0-10% 


748 ± 093 


i-.OOiJ u.u±± 


— n ni 8 

u.uxo 


Au+Au 


70-80% 


946 ± 140 


0.980 + 0.015 


- 0.027 


Au+Au 


60-70% 


964 ± 081 


0.9930 + 0.0036 - 


- 0.0059 


Au+Au 


50-60% 


0.971 + 0.053 


0.9965 + 0.0011 - 


- 0.0019 


Au+Au 


40-50% 


0.959 ± 0.053 


0.9982 + 0.0008 - 


- 0.0012 


Au+Au 


30-40% 


0.956 + 0.051 


0.9991 + 0.0005 - 


- 0.0008 


Au+Au 


20-30% 


0.969 ± 0.035 


0.9996 + 0.0004 - 


- 0.0006 


Au+Au 


10-20% 


0.980 + 0.021 


0.99998 + 0.00019 - 


0.00026 


Au+Au 


5-10% 


0.976 + 0.020 


1.00003 + 0.00011 - 


0.00019 


Au+Au 


0-5% 


0.986 + 0.013 


1.00023 + 0.00008- 


0.00012 


Au+Au ZDC 


"20-30%" 


0.944 ± 0.057 


0.99940 + 0.00070 - 


0.00084 


Au+Au ZDC 


"10-20%" 


0.980 ± 0.022 


0.99963 + 0.00010 - 


0.00012 


Au+Au ZDC 


"5-10%" 


0.982 + 0.015 


1.00001 + 0.00005 - 


0.00009 


Au+Au ZDC 


0-5% 


0.994 + 0.005 


1.00021 + 0.00003 - 


0.00005 
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determined by the 3-particle ZYAM, using 10% of the bins, with the default value 
of b. We then change the number of bins used in the ZYAM from 10% to 5% and 
15%. Figure 13.271 shows the background subtracted 3-particle correlations when b is 
obtained from the 3-particle ZYAM using the default a for 5% and 15% of the bins. 
Table [37^ shows the default values for b and the systematic uncetainties applied. The 
uncertainties listed are only from this change in the number of bins. The change in b 
that is fully correlated with the uncertainty in a is not listed (i.e. how much b changes 
to perseve 3-particle ZYAM for variations of a within its uncertainty). 

3.5.2 Flow Correlation 

There are several systematics associated with flow correlations. The dominant 
systematic from flow is the uncertainty on the flow measurement. The other sources 
of systematic error due to flow are much smaller. These include systematics associated 
with using the trigger particle flow for the jet flow, effect of V2 fluctuations, and the 
trigger particle v^. 

There are several measurements of elliptic flow: the reaction plane method, the 
modified reaction plane method, the 2-particle cumulant and the 4-particle cumulant. 
The different measurements give different results because the measured V2 have dif- 
ferent sensitivity to other effect such as non-flow. The reaction plane method over 
estimates the V2 due to non-flow. The modifled reaction plane method has reduced 
sensitivity from non-flow at high pt, and will not over estimate the flow as much 
as the standard reaction plane method. The 2-particle V2 over estimates the flow 
due to non-flow and f 2 fluctuations. The 4-particle cumulant has little sensitivity to 
non-flow and will underestimate the flow in the presence of V2 fluctuations. For our 
default V2 value we use the average of the modifled reaction plane and the 4-particle 
measurement. For our systematics we use the 4-particle cumulant V2 for the lower V2 
value and the 2-particle cumulant t>2 value is used for the upper limit. The 2-particle 
cumulant is used as the upper limit instead of the modifled reaction plane because 



Figure 3.27. Background subtracted 3-particle correlations for the un- 
certainty on normalization factor b. Left: b from 5% of bins for ZYAM. 
Center: b from 15% of the bins for ZYAM. Right: Off-diagonal away-side 
projection from default with systematics uncertainty from the uncertainty 
on b shown in shaded band. From top to bottom plots are pp, d+Au, 
Au+Au 50-80%, Au+Au 30-50%, Au+Au 10-30%, Au+Au 0-10%, and 
ZDC triggered Au+Au 0-12% coUisions at = 200 GeV/c. 
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the 2-particle cumulant includes the systematic uncertainty due to V2 fluctuations. 
In the most central and the most peripheral bins we do not have measurements for 
the 4-particle cumulant V2- The 4-particle cumulant V2 is estimated from extrapola- 
tion to be about 50% of the reaction plane V2 (see Fig. 12.21) . Figure [3^281 shows the 
background subtracted 3-particle correlations using the 2-particle cumulant V2 and 
4-particle cumulant V2. 

One of the other flow systematics is on the jet-fiow. By default the trigger- 
associated pair is assumed to have the same azimuthal anisotropy as the trigger 
particle. To assess the effect of this, we assume that the azimuthal anisotropy of the 
trigger-associated pair may be more like the azimuthal anisotropy of a particle with 

use the V2 of 5 GeV/c particles which is approximately the total 
Pt of a trigger particle of 3 < < 4 GeV/c and an associated particle of 1 < p-r < 2 
GeV/c to estimate the systematic uncertainty due to jet-fiow. The jet-fiow is put into 
the analysis in the construction of the hard-soft term by folding the flow modulated 
background with jet-fiow value. The background subtracted 3-particle correlation 
using the 5 GeV/c particle flow is shown in Fig. 13.291 along with the default jet-fiow 
for comparison. As an additional cross check, the background subtracted 3-particle 
correlation was examined without the jet-fiow subtraction. This was done as a check 
to see if observing a conical emission signal was dependent on whether on not this 
subtraction was performed. It is not included in the systematic error bars. The signal 
is reduced but still present without this subtraction as shown in Fig. 13.291 right. 

The estimate of this systematic uncertainty on the azimuthal anisotropy of the 
trigger associated pair can be improved in the future by using 2-particle correlations 
for trigger particles with different orientations with respect to the reaction plane. 
The 2-particle correlation can be divided into bins depending on the trigger particle 
orientation with respect to the reaction plane. This would provide the reaction plane 
dependence of the trigger associated pair. The hard-soft term could be constructed 
for each of these bins and summed together with a number of trigger particle weight- 
ing. This would provide a better estimate of the azimuthal anisotropy of the trigger 



Figure 3.28. Background subtracted 3-particle correlations for uncertainty 
on V2- Left: V2 from the 2-particle cumulant measurement. Center: V2 
from the 4-particle cumulant method. Right: Off-diagonal away-side pro- 
jection from default with systematic uncertainty from the uncertainty on 
V2 shown in shaded band. From top to bottom plots are Au-|-Au 50- 
80%, Au+Au 30-50%, Au+Au 10-30%, Au+Au 0-10%, and ZDC triggered 
Au-FAu 0-12% collisions at ^/s]^ = 200 GeV/c. 
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associated pair than what is currently used. Once this has been performed, since it is 
a better estimate for the jet-flow, this should probably be used for the default value 
and the current value should be used for the assessment of the systematic uncertainty. 

Another uncertainty on the flow is from the parameterization of the ^4. The V4 
was fit to the ratio of the f4/f| for 1 < p-p < 2 GeV/c. This was not done for our 
trigger particles of 3 < < 4 GeV/c because the f4 data does not go out to this 
high in pt- To account for the uncertainty on the trigger particle W4 due to this a 
±20% uncertianity on the trigger particle f 4 has been applied. Figure 13.301 shows the 
background subtracted 3-particle correlations with the trigger particle V4 decreased 
by 20% and increased by 20%. 

There is an additional uncertainty on the elliptic flow in the ZDC triggered data. 
There have been no flow measurements made for the ZDC data so we have used the 
minimum bias flow values. These have been applied such that a ZDC event uses the 
V2 value from the centrality bin that it would fall into if it were a minimum bias 
event. A cross check has been applied by parameterizing the V2 values as a function 
of multiplicity and using the value from the parameterization. Due to first order 
cancellation of the V2 * f 2 and * V4 this was found to have a negligible impact on 
the final result. 

3.5.3 Finite Multiplicity Bin Systematics 

There can be effects on the soft-soft and hard-soft background terms from doing 
the analysis in finite multiplicity bins. In the soft-soft background term, the pedestal 
and the A0i, A02 dependent structure scale differently with multiplicity. As an 
assessment of the uncertainty in the soft-soft term from using a finite width centrality 
bin due to this difference in scaling, the soft-soft term was constructed in individual 
multiplicity bins. The individual soft-soft terms were summed together in a weighted 
sum to make the soft-soft term for the entire multiplicity bin. They were weighted 
by an estimated underlying event multiplicity distribution. The underlying event 
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Figure 3.29. Background subtracted 3-particle correlations for uncertainty 
on the jet-flow. Left: Jet-flow from 5 GeV/c particles. Left Center: Jet- 
flow same as trigger particle, default. Right Center: Off-diagonal away- 
side projection from default with systematic uncertainty from the uncer- 
tainty on jet flow shown in shaded band. Right: No jet-flow subtraction. 
From top to bottom plots are Au+Au 50-80%, Au-^Au 30-50%, Au-^Au 
10-30%, Au-fAu 0-10%, and ZDC triggered Au+Au 0-12% collisions at 
= 200 GeV/c. 
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Figure 3.30. Background subtracted 3-particle correlations for uncertainty 
on the trigger particle f4. Left: Trigger particle ^4 decreased by 20%. Cen- 
ter: Trigger particle f4 increased by 20%. Right: Off-diagonal away-side 
projection from default with systematic uncertainty from the uncertainty 
on trigger particle V4, shown in shaded band. From top to bottom plots 
are Au+Au 50-80%, Au+Au 30-50%, Au+Au 10-30%, Au+Au 0-10%, 
and ZDC triggered Au-^Au 0-12% collisions at y/s^ = 200 GeV/c. 
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multiplicity distribution was estimated from the trigger event multiplicity distribution 
minus the number of jet-like correlated particles counted in the multiplicity. The 
number of jet-like correlated particles was estimated from a background subtracted 
2-particle correlation using trigger particles with the same cuts as in the analysis and 
associated particles with the reference multiplicity cuts. However, due to the way 
we have stored the subset of the data we process the number of fit points cut was 
15 instead of 10. This estimate had to be rounded to the nearest integer since the 
multiplicity distributions are integer. Figure [3.311 shows the the effect of using the soft- 
soft term constructed this way on the background subtracted 3-particle correlation 
results. 

The hard-soft term can also have effects due to the finite multiplicity bin width. 
Optimally the term would be created on an event-by-event basis; however, this cannot 
be done. To check the magnitude of the effect, the ZDC trigger Au+Au data was 
dividing into multiplicity bins 1/10*^ the size. This effect should be the largest in 
central collisions where the largest change in the hard-soft term was 5% that of the 
peaks, as shown in Fig. 13.321 This could not be done in bins I/IO*'' the size of the 
default bins for all centralities and data sets due to statistics, so an overall ±5% 
systematic uncertianity has been applied. This systematic can be improved in the 
future by dividing the each of the centrality bins into as many subbins as the statistics 
will allow. This will take into account the structure of this correction on A(/)i, A02 
space. 

3.6 Summary 

Three-particle azimuthal correlations have been measured to discover the physics 
mechanism(s) behind the broadened and maybe even double-peaked away-side struc- 
ture in 2-particle correlations in central Au+Au collisions. Three-particle correlations 
have the power to distinguish physics mechanisms with conical emission (hydrody- 
namic conical flow generated by Mach-cone shock waves and Cerenkov gluon radia- 
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Figure 3.31. Background subtracted 3-particle correlations for uncertainty 
on the soft-soft term. Left: Soft-soft term constructed from weighed sum 
of soft-soft terms for individual multiplicities. Center: Default. Right: 
Off-diagonal away-side projection from default with systematic uncer- 
tainty from the uncertainty on soft-soft shown in shaded band. From 
top to bottom plots are Au+Au 50-80%, Au-hAu 30-50%, Au-hAu 10- 
30%, Au+Au 0-10%, and ZDC triggered Au+Au 0-12% collisions at 
= 200 GeV/c. 
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Figure 3.32. Background subtracted 3-particle correlations for uncertian- 
ity on the hard-soft on the hard-soft term from bin width effects in ZDC 
triggered Au-|-Au 0-12%. Top Left: Default. Top Center: With hard-soft 
from finer multiplicity bins. Top Right: Difference. Bottom: On-diagonal, 
left, and off-diagonal projections, right. The points are from the default 
and the histograms are from the finer mulitplicity binned hard-soft term. 



tion) from other physics mechanisms. This analysis is designed to extract the jet-like 
3-particle correlations, by treating the event as composed of two components, particles 
jet-like correlated with the trigger particle and particles not jet-like correlated with 
the trigger particle (but correlated via anisotropic flow). We assume the background 
subtracted 3-particle correlations are positive deflnite by normalizing the background 
via 3-particle ZYAM. The systematics on this analysis have been studied in great 
detail. 

The 3-particle azimuthal correlations have been studied for 'p'p^ d+Au, and dif- 
ferent centralities of Au+Au collisions (for both minimum bias and ZDC triggered 
central data) for a trigger particle of 3 < p^^*^ < 4 GeV/c with two associated particles 
of 1 < < 2 GeV/c. We observed signiflcant off-diagonal peaks, the signature 

of conical emission, in the mid-central and central Au-|-Au collisions at about 1.42 
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radians from tt. In ZDC triggered central Au+Au collisions, the associated particle 
Pt dependence has been studied. No strong pt dependence of the emission angle has 
been observed, suggesting the observed conical emission is due to Mach-cone shock 
waves, not Cerenkov gluon radiation. The average speed of sound of the medium is 
extracted to be Cg = 0.15c in a simple-minded Mach-cone scenario. Further and more 
refined theoretical studies are urgently needed to assess the effects of hydrodynamic 
expansion and different scenarios of QGP-hadron phase transition and to connect our 
measurement to the properties of the medium and its equation of state. 
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4. CONCLUSIONS AND RECOMMENDATIONS 

Jet quenching-suppression of high yields and correlations-observed at RHIC con- 
stitutes important evidence for the discovery of the created medium being a nearly 
perfect fluid. In particular, jet-like correlations both at high pt and low pt have 
provided a valuable tool to study the properties of the medium. This thesis con- 
tributes to the enterprise of jet-like correlation studies at RHIC. Particularly, the 
3-particle correlations, presented in this thesis, are the first such measurement and 
should provide new and unique insights into the medium at RHIC 

4.1 Two-Particle Correlations of Charged Hadrons 

Jet-like correlations of charged hadrons have been studied in d+Au collisions. 
This is the first d+Au jet-like correlation analysis to go down to very low associated 
Pt- The correlations were found to be consistent with previous pp results. This 
demonstrates that it is the final state in Au+Au collisions, not the initial state of 
Au nuclei, that is responsible for the modifications to jet-like correlations observed 
in central Au+Au collisions. This also lets us use the d+Au data besides pp as a 
reference for the Au+Au data to gain the benefit of the increased statistics in d+Au. 
The widths of the correlation functions and the spectra of the correlated yields have 
been analyzed for d+Au collisions and can be used for comparison to theoretical 
calculations. 

Jet-like correlations of charged hadrons have also been studied in Au+Au col- 
lisions. This was done using Au+Au data from the fourth year of RHIC running 
which gave an increase in statistics of about an order of magnitude from the second 
year data. This analysis was done using the same trigger and associated pt ranges 
1.0 < < 2.5 < Pt"^ < 4 GeV/c as for the published PHENIX resuhs [48J. 
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The PHENIX data show a significant dip at tt in the 2-particle azimuthal correlations 
that is strongest in mid-central Au+Au coUisions. Our results also show a dip at tt; 

however, the dip in our results is not as strong. In our results, the dip is strongest in 
the most central collisions, but is not very significant given the systematic uncertainty 
on the elliptic flow measurement. 

4.2 Two-Pcirticle Correlations with Identified Trigger Particles 

Jet-like correlations with identified trigger particles have been studied in d-l-Au 
collisions and Au+Au collisions in difi^erent centralities. Proton-proton collisions did 
not have adequate statistics for this study. The interest in studying identified cor- 
relation is due to the baryon/meson puzzle, the large baryon over meson ratio at 
intermediate px in central Au-|-Au collisions. The analysis was done for intermediate 
Pt trigger particles where the large baryon/meson ratio is observed and where we 
have good statistics and particle identification capability. Trigger particles are iden- 
tified as p, p, TT"*", and 7r~ by the relativistic rise of dE/dx. No significant difference is 
observed between correlations with identified baryon and meson triggers. The large 
baryon/meson ratio does not seem to have a large effect on the particles correlated 
with the baryons and mesons. This imposes a serious challenge to the initial coa- 
lescence and recombination models which successfully and elegantly explained the 
large baryon/meson ratio at intermediate pt and the constituent quark scaling of the 
elliptic flow. 

For this analysis to be flnalized, a detailed study of the systematics needs to be 
carried out. Since the V2 correlation will partially cancel in the difference between 
baryon and meson triggered correlations, the systematic effects of flow subtraction 
will not be very large. The systematic effects due to background normalization should 
cancel to flrst order and be negligible in the difference. A study needs to be done 
on the purities. The purities depend on the measured to charged tt ratio. This 
has two systematics associated with it. The flrst is the assumption that the yields 
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for — K~ — Kg. However, this systematic effect should be neghgible if one sums 
p+p and 7r+ + 7r~ as K'^ + K~ — 2K'^ to a good approximation. The second is the 
systematics due to the error in the A'^ measurement. 

Other similar analyses can be done that are of interest. One is to use trigger 
particles identified through VO reconstruction such as A, and K^. This analysis is 
already on-going in STAR. Another is to use a higher pt hadron trigger and use 
the identified particles as our associated particles. This would allow one to look for 
baryon and meson dependences on the associated particles. It would also be beneficial 
to look at both identified trigger and associated particles, at both low and high px, 
to see if there are differences in baryon-baryon, baryon-meson, meson-baryon, and 
meson-meson correlations. Using charged hadrons as the trigger and identifying the 
associated particles may be an analysis that can be done with the current data. Both 
analyses (and the current analysis) would benefit from increased statistics. PHENIX 
has done a good job in identified particle correlations using their time of fiight detector 
for particle identification. Their results are currently limited to intermediate pr- The 
Pt dependences of the identified correlations will be important to shed light on some 
of the puzzles. There will be improvements to particle identification in STAR with 
the inclusion of the time of fiight detector that will greatly benefit all of the identified 
particle analyses. 

4.3 Three-Pcirticle Correlations 

Jet-like 3-particle azimuthal correlations have been studied in pp, d-l-Au, and dif- 
ferent centralities of Au-|-Au collisions. This analysis was done to investigate the 
physics mechanism(s) behind the broadened or double-peaked away-side correlation 
structure observed in central Au-|-Au collisions. These physics mechanisms include: 
large angle gluon radiation, jet deflection by radial flow or preferential selection of 
particles due to path-length dependent energy loss, hydrodynamic conical flow gener- 
ated by Mach-cone shock waves, and QCD Cerenkov gluon radiation. Three-particle 
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correlations can identify conical emission which can be generated by Mach-cone shock 
waves or Cerenkov gluon radiation, but not by other proposed physics mechanisms. 
The analysis was carried out treating the triggered event as composed of particles 
jet-like correlated with the trigger and background particles to extract the genuine 
3-particle jet-like correlations. The combinatorial background was normalized assum- 
ing the 3-particle jet-hke correlation is positive definite using ZYAM (zero yield at 
minimum). The systematics have been studied in great detail. 

We have observed significant off-diagonal peaks-the signature for conical emission- 
in mid-central and central Au-|-Au collisions. These peaks are found to be at about 
TT ± 1.42 radians from the trigger particle (i.e. emission angle). We have studied the 
associated particle pr dependence and found the emission angle to have no significant 
dependence on the associated particle pr- This suggests that the observed conical 
emission is due to Mach-cone shock waves and not due to Cerenkov gluon radiation, 
which predicts an angle that is sharply decreasing with associated particle pr- 

Several recommendations can be made for future studies of 3-particle correlations. 
They include: 

• Trigger particle px dependence. This thesis has attempted this study but could 
not carry it to very high px due to statistics. By increasing the trigger particle 
Pt one can probe different relative strengths of conical emission and back-to- 
back emission. Also correlations with higher pt trigger benefit from an increased 
signal to noise ratio due to an increase in the signal. Given enough statistics, 
potentially from an online high px or high trigger, the increased signal to 
noise ratio could greatly reduce the systematic errors. 

• Three-particle correlations with identified associated particles (such as p and 
tt). This can provide an important check on Mach-cone emission since protons 
and pions are expected to have different associated particle px dependence on 
the emission strength. 
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• Three-particle correlations with two trigger particles and one associated particle 
to probe different jet samples by varying the pt of the two triggers. 

• Trigger particle species dependence to look for differences between hght quark 
triggers, heavy quark triggers, and non-photonic electron triggers. 

• The 3-particle correlation study in this thesis has concentrated on the away 
side. There is also an interesting phenomenon on the near side. The near-side 
ridge-particles that are jet-like correlated with the trigger particle in but not 
in r^-can be studied through 3-particle Ary-Ar/ correlations. 

In addition, 4-particle azimuthal correlations can be explored. This could be done 
using two trigger particles and two associated particles. The second trigger particle 
would allow for the determination of the 77 of the away-side jet-axis (which is not 
correlated with the near-side). The two associated particles can then be studied in 
a polar coordinate system determined by the second trigger particle. This is a more 
natural coordinate system to study conical emission. Such an analysis would be very 
complicated due to a large number of background terms that can each be rather 
complicated. 
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APPENDIX 

Figure lA.ll shows the 0-averaged detector efficiency for 3 centrahties of d+Au coUi- 
sions and 9 centrahties of Au+Au cohisions. The efficienies are fit to, 

pOe ^PT> . (A.l) 

The inverse of the fit is used to correct for the number of particles at the single particle 
level. 
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Figure A.l. Detector efficiency for charged pions for 3 centralities of d+Au 
collisions (top left) and 9 centralities of Au+Au collisions as a function of 
Pt at ^/sNN = 200 GeV. 



109 



Figures IA.2I and IA.3I show Na.,, distributions in Au+Au and d+Au respectively. 
The distributions are fit to 6 Gaussians, one each for p, vr"*", vr", , K~ . The 
yields for tt^, tt^, p, and p are pO, p2, p3, p5, respectively. The positive and negative 
particles share the same centroid. The centroids are given by pi minus 10 and p4 
minus 10 for ti^ and p/p, respectively. All peaks share a common width given by p7. 
The yields are fixed by the measured to tt ratio. The are fixed to 2/3 of 
the distance between the pion and proton centroid (closer to the proton's centroid) 
which was determined from the Bethe-Bloch formula. 



Figure A. 2. Na-,, distributions with fits in Au+Au collisions at ^Jsnn = 200 
GeV/c. Curves are from fit with red for pion, blue for proton and green 
for kaon. Rows are the centrality bins 70-80%, 60-70%, 50-60%, 40-50%, 
30-40%, 20-30%, 10-20%, 5-10% and 0-5% from top to bottom. Columns 
are for different pt and rj cuts. Left: |?7| < 0.35 and 3 < pt < 3.5 GeV/c. 
Left center: 0.35 < \rj\ < 0.75 and 3 < pt < 3.5 GeV/c. Right center: 
\ri\ < 0.35 and 3.5 < pt < 4 GeV/c. Right: 0.35 < |r/| < 0.75 and 
3.5 < pt < 4: GeV/c. Fit parameters are p0=7r~ yield, pl=7r centroid-|-10, 
p2=7r+ yield, p3=p yield, pA=p/p centroid+10, p5=p yield, and p7=width 
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Figure A. 3. Same as Fig. 6.2 but for d+Au. Rows are the centrality bins 
100-40%, 20-40% and 0-20% from top to bottom. 



Figures IA.4tiA.7l show the Not^ distributions for positive and negative particles 
seperately with hne representing cuts. Particles with A^cr,r greater than the right line 
are vr^ with a 95% purity. Particles with Na-,^ less than the left line are with a 
50% purity. The curves are from the fits shown in Figures [A. 21 and IA.3I 
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Figure A. 4. Na-j^ distributions with fits and cuts for positive particles in 
Au+Au collisions at yj/vw = 200 GeV/c. Curves are from fit with red for 
pion, blue for proton and green for kaon. Particles to the left of the left 
line are 50% protons and particles to the right of the right line are 95% 
pion. Rows are the centrality bins 70-80%, 60-70%, 50-60%, 40-50%, 30- 
40%, 20-30%, 10-20%, 5-10% and 0-5% from top to bottom. Columns are 
for different px and rj cuts. Left: |?7| < 0.35 and 3 < pr < 3.5 GeV/c. Left 
center: 0.35 < \t]\ < 0.75 and 3 < pr < 3.5 GeV/c. Right center: |?7| < 0.35 
and 3.5 < px < 4.0 GeV/c. Right: 0.35 < \ri\ < 0.75 and 3.5 < pr < 4.0 
GeV/c. 
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Figure A. 5. Same as Fig. 6.4 but for negative particles. 
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Figure A. 6. Same as Fig 6.4 but for d+Au. Rows are the centrality bins 
40-100%, 20-40% and 0-20% from top to bottom. 




NSIgma Plan HSIgma Plan HSigiiaPIgn NSIgma Plan 



Figure A. 7. Same as Fig. 6.6 but for negative particles. 

The A^o"^ cuts on the trigger particles are fit as a function of trigger particle pt 
to reduce the error on the fit. Figures [A. SUA. 101 show these fits. Although we fit the 
A^cTtt out to higher px only trigger particles of 3 < p^r < 4 GeV/c were used, due to 
limited statistics at high px- The cuts are fit to a second order polynomial. 
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Figure A. 8. Fits to the cuts on Na-^ as a function of pt for Au+Au collisions 
at y/SNN = 200 GeV/c. Rows correspond to the centrality bins 70-80%, 
60-70%, 50-60%, 40-50%, 30-40%, 20-30%, 10-20%, 5-10% and 0-5% from 
top to bottom. Columns are for differents charges and 1] cuts of pions. Left: 
vr"*" with \ri\ < 0.35. Left center: tt"*" with 0.35 < \ri\ < 0.7. Right center: 
TT" with \ri\ < 0.35. Right: vr" with 0.35 < \r]\ < 0.7. 
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Figure A. 9. Same as Fig. 6.8 but for p{p). Left: p with |?7| < 0.35. Left 
center: p with 0.35 < |?7| < 0.7. Right center: p with |?7| < 0.35. Right: p 
with 0.35 < \r]\ < 0.7. 
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Figure A. 10. Same as Fig. 6.8 but for d+Au. Rows correspond to the 
centrality bins 40-100%, 20-40%, and 0-20% from top to bottom. 
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Figure A. 11. Same as Fig. 6.9 but for d+Au. Rows correspond to the 
centrahty bins 40-100%, 20-40% and 0-20% from top to bottom. 



Figures IA.12tlA.19] show the 0-dependence of the acceptance. This ^-dependence 
is due to the TPC sector boundries. The plots are normahzed such that the average 
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is one since the 0-average efficiency is obtianed seperately (Fig. lA.ip . Trigger and 
associated particles were corrected for this 0-dependence on the single particle level. 




Figure A. 12. TPC acceptance in for particles of 1 < < 2 GeV/c 
with +0.5 Tesla magnetic field. Positive particles are shown in black and 
negative particles are shown in red. Plots (left to right and top to bottom) 
are for Au+Au collisions at a/s/vtv = 200 GeV/c in centralities 70-80%, 
60-70%, 50-60%, 40-50%, 30-40%, 20-30%, 10-20%, 5-10%, and 0-5%. 
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Figure A. 13. Same as Figure 6.12 but for —0.5 T magnetic field. 




Figure A. 14. Same as Fig. 6.12 but for 3 < < 4 GeV/c. 




Figure A. 15. Same as Fig. 6.12 but for 3 < < 4 GeV/c and —0.5 T 
magnetic field. 
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Figure A. 16. Same as Fig. 6.12 but for ZDC trigger central Au+Au col- 
lisions in centralities (from left to right and top to bottom) "20-30%", 
"10-20%", "5-10%", and 0-5% most central where the first three centralites 
do not correspond to the actual percantage of cross section due to the ZDC 
trigger but have the same multiplicty cuts as the minimum bias events that 
do correspond to the given cross section. 




Figure A. 17. Same as Fig. 6.16 but for —0.5 T magnetic field. 
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Figure A. 18. Same as Fig. 6.16 but for 3 < < 4 GeV/c. 
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Figure A. 19. Same as Fig. 6.16 but for 3 < < 4 GeV/c and —0.5 T 
magnetic field. 
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Figure A. 20. Explanation of panels for 3-particle analysis plots. 



Figure IA.20I is a graphical explanation of the various panels in Figures IA.2HlAl34l 
These figures are the working plots for the 3-particle analysis in Chapt. 4. In these 
figures, the top left and top left center panels show the raw 3-particle correlation 
function in two different representations. The top center and top right center show the 
hard-soft background term. The top right panel shows the raw 2-particle correlation 
function in red, the mixed event background with flow modulation from V2 and V4 
in solid black, and the normalized to 3-particle ZYAM mixed event background with 
flow modulation in open black. In the second row, the left and left center plots shown 
the raw 3-particle correlation function with the hard-soft background subtracted. 
The center and right center plots show the soft-soft background term. The right 
plots show the background subtracted 2-particle correlation function in red, the flow 
from V2 in black, and the flow from V4 in blue. In the third row, the left and left 
center plots show the raw 3-particle correlation function with the soft-soft background 
subtracted. The center and right center panel shows the flow contribution from elliptic 
flow, V2, between the trigger and associated particles. The right panel shows the flow 
background from the non-flow structure in the soft-soft term flowing with the trigger 
particle. In the fourth row, the left and left center panels show the additional 3- 
particle flow between the trigger and associated particle when contributions are 
considered. The center and right center panels show the normalized sum of the 
soft-soft and flow terms. The right panel is an off-diagonal projection of the soft-soft 
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term. In the bottom row, the left and left center panels show the sum of the hard-soft 
terms and the flow terms. The center and right center panels show the background 
subtracted 3-particle correlation. The bottom right plot is not used. 




Figure A. 21. Three-particle correlation analysis plots for 3 < p^^*^ < 4 
GeV/c and 1 < < 2 GeV/c in pp collisions at ^/snn = 200 GeV/c. 
Explanation of panels is in preceeding table and text. 




Figure A. 22. 



Same as Fig. 6.21 but for d+Au collisions. 
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Figure A. 23. Same as Fig. 6.21 but for 50-80% Au+Au collisions. 
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Figure A. 24. Same as Fig. 



6.21 but for 30-50% Au+Au collisions. 
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Figure A. 25. Same as Fig. 6.21 but for 10-30% Au+Au collisions. 
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Figure A. 26. 



Same as Fig. 6.21 but for 0-10% Au+Au collisions. 
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Figure A. 27. Same as Fig. 6.21 but for 0-12% ZDC triggered Au+Au 
collisions. 
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Figure A.29. Same as Fig. 6.27 but for 0.75 < < 1.0. 
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Figure A.30. Same as Fig. 6.27 but for 1.0 < pj^''"" < 1.5. 
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Figure A.31. Same as Fig. 6.27 but for 1.5 < pj^''"" < 2.0. 
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Figure A.32. Same as Fig. 6.27 but for 2 < pj,''"'' < 3. 
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Figure A. 33. Same as Fig. 6.27 but for 4 < p^"^ < 6. 
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Figure A. 34. Same as Fig. 6.27 but for 4 < p^"^ < 6. 

Figure lA. 351 shows the background subracted 3-particle correlations in the central- 
ity bins in which they are analyzed. The bins used for the minimazition are outlined 
in black. These bins are summed together with a number of trigger particle weight- 
ing to consturct the background subtracted 3-particle correlation plots in the wider 
centrality bins. 



139 




Figure A. 35. Background subtracted 3-particle correlations with the bins 
used for the 3-particle ZYAM highlighted in black. Panels are, from left 
to right and top to bottom, pp 0-100%, d+Au 20-100%, 10-20%, and 0- 
10%, Au+Au 70-80%, 60-70%, 50-60%, 40-50%, 30-40%, 20-30%, 10-20%, 
5-10%, and 0-5% and central ZDC triggered Au+Au "20-30%", "10-20%", 
"5-10%", and 0-5% coUisions at = 200 GeV/c. Where the cross 

sections in quotes are not the actual cross sections but the minimum bias 
cross section cuts on the centeral triggered data. 



Figrue IA.36I shows the on-diagonal and off-diagonal projections with the Gaussian 
fits. Each projection is fit to a central Gaussian and symmetric side Gaussians. The 
parameters listed for the fits are (from top to bottom) x^/ndf, side Gaussian yield, 
side Gaussian distance from the center (in radians), side Gaussian width (in radians), 
central Gaussian yield, and central Gaussian width. The central Gaussian is centered 
at IT (on-diagonal) or zero (off-diagonal). 
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Figure A. 36. Backgroud subtracted 3-particle correlations shown with on- 
diagonal (center) and off-diagonal (right) projections of the away-side dis- 
tribution. The curve represent fits to a central Gaussian and symmetric 
side Gaussians. The rows are (from top to bottom) pp, d-l-Au, Au+Au 50- 
80%, Au+Au 30-50%, Au+Au 10-30%, Au+Au 0-10% and ZDC triggered 
Au+Au 0-12% colhsions at ^/snn = 200 GeV/c. The projections are of a 
strips of full width of 0.7 radians. Errors are statistical. 
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Figures IA.37I and IA.38I show the on-diagonal projections of the away-side of the 
background subtracted 3-particle correlations. The on-diagonal projections are shown 
here with systematic errors. 




(A(|)+A(|)^)/2 (A*/A^^)/2 (A(|)+A(l)^)/2 



Figure A. 37. On-diagonal projections of the background subtracted 3- 
particle correlations in strips of full width 0.7 radians on the away-side. 
From left to right, top to bottom are pp, d+Au, Au+Au 50-80%, Au+Au 
30-50%, Au+Au 10-30%, Au+Au 0-10% and ZDC triggered Au+Au 0-12% 
collisions at ^/snn = 200 GeV/c. Systematic errors are shown in the yellow 
boxes. Dashed black lines are at zero. 



Figures IA.39I and IA.40I show the on-diagonal and off-diagonal projections of the 
near-side of the background subtracted 3-particle correlations. The projections are 
shown here with systematic errors. 
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Figure A. 38. Same as previous exept panels are from left to right and 
top to bottom, 0.5 < p^''°^ < 0.75 GeV/c, 0.75 < < 1.0 GeV/c, 

1.0 < p^''°^ < 1.5 GeV/c, 1.5 < p^"""" < 2.0 GeV/c, and 2.0 < p^''°^ < 3.0 
GeV/c for 0-12% ZDC triggered Au+Au. 
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Figure A. 39. On-diagonal projections (blue) and off-diagonal projection 
(red) of the background subtracted 3-particle correlations in strips of full 
width 0.7 radians on the away-side. From left to right, top to bottom are 
pp, d+Au, Au+Au 50-80%, Au+Au 30-50%, Au+Au 10-30%, Au+Au 0- 
10% and ZDC triggered Au+Au 0-12% collisions at = 200 GeV/c. 

Systematic errors are shown with the blue histograms and yellow boxes for 
the on-diagaonl and off-diagonal projections, respectively. Dashed black 
lines are at zero. 



Figure IA.41I shows the on-diagonal and off-diagonal projections of the near-side 
peak. The lines show fits to a single Gaussian centered at zero. The parameters listed 
for the fits are (from top to bottom) x^/ndf, Gaussian yield, and Gaussian width (in 
radians) . 
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Figure A. 40. Same as previous exept panels are from left to right and 
top to bottom, 0.5 < vf'""^ < 0.75 GeV/c, 0.75 < < 1.0 GeV/c, 

1.0 < p^''°^ < 1.5 GeV/c, 1.5 < p^"""" < 2.0 GeV/c, and 2.0 < p^''°^ < 3.0 
GeV/c for 0-12% ZDC triggered Au+Au. 



145 



|3.0<PI<4.01.0Qcpt<2.00100-PEl | | lxty),'2 with Wldth^O.7 




Figure A. 41. Backgroud subtracted 3-particle correlations shown with on- 
diagonal (center) and off-diagonal (right) projections of the near-side peak. 
The curve represent fits to a Gaussian centered at zero. The rows are 
(from top to bottom) pp, d+Au, Au+Au 50-80%, Au+Au 30-50%, Au-^Au 
10-30%, Au+Au 0-10% and ZDC triggered Au+Au 0-12% coUisions at 
■s/snn = 200 GeV/c. The projections are of a strips of full width of 0.7 
radians. Errors are statistical. 
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